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Abstract 

The  deleterious  effects  of  fretting  on  the  fatigue  properties  of  a  material  have  been  known  since 
the  work  of  Warlow-Davies  in  1941.  However,  a  widely  accepted  fretting  fatigue  life  prediction 
method  still  does  not  exist  and  debates  persist  as  to  what  the  critical  parameter(s)  for  fretting 
fatigue  are.  This  work  demonstrates  that  the  surface  stresses  due  to  contact  can  be  used  to 
characterize  the  damaging  effects  of  fretting  on  the  fatigue  properties  of  a  material. 

The  elastic  stress  analysis  performed  by  Hamilton  for  a  sphere  on  flat  plane  fretting  geometry 
is  used  to  determine  the  stresses  due  to  contact  for  various  experimental  conditions  provided 
in  previous  work.  It  is  shown  that  the  specimen  fretting  fatigue  life  is  directly  related  to  the 
surface  stresses  induced  in  the  material  by  fretting. 

Fretting  fatigue  life  prediction  methods  are  analyzed.  The  stress-life  approach,  which  is 
widely  used  in  industry,  is  examined  for  its  applicability  for  fretting  fatigue  conditions.  The 
damage  tolerant  approach  is  used  to  determine  the  relative  magnitudes  of  the  crack  nucleation 
and  long  crack  propagation  stages  of  fretting  fatigue  life. 

An  analysis  of  fretting  fatigue  palliatives,  focusing  on  shot  peening  and  coatings,  is  pre¬ 
sented.  A  systematic  method  to  determine  the  optimum  shot  peening  depth  for  fretted  com¬ 
ponents  is  proposed.  The  method  uses  the  elastic  stress  field  expressions  derived  by  Hamilton 
for  a  sphere  on  flat  surface  contact  geometry  and  is  demonstrated  for  Ti-6A1-4V.  However,  the 
method  is  general  and  can  be  applied  for  any  material  and  any  contact  geometry  which  allows 
analytic  evaluation  of  the  stress  fields. 

Thesis  Supervisor:  Subra  Suresh 

Title:  Richard  P.  Simmons  Professor  of  Materials  Science  and  Engineering 
Professor  of  Mechanical  Engineering 
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Chapter  1 


Introduction 


1.1  Definitions  and  Background 

Fretting  fatigue  generally  refers  to  the  degradation  of  the  fatigue  properties  of  a  material  due 
to  cyclic  relative  displacement  between  two  contacting  materials.  The  most  severe  impact  on 
the  fatigue  properties  occurs  when  displacement  amplitudes  are  10-50  fim  [21].  The  fretting 
action  produces  considerable  wear  damage  to  the  material  and  also  produces  very  large  contact 
stresses.  The  reduction  in  fatigue  life  under  these  circumstances  can  be  quite  significant.  The 
National  Materials  Advisory  Board  [l]  gives  the  following  definitions. 

•  Fretting  .  Fretting  is  a  surface  damage  phenomenon  occurring  on  two  contacting  surfaces 
having  oscillatory  relative  motion  of  small  amplitude. 

•  Fretting-Initiated  Fatigue  .  A  condition  where  the  material  fatigue  strength  is  de¬ 
graded  by  the  presence  of  the  following: 

1.  Surface  stress  concentrations  resulting  at  sites  of  fretting  pits. 

2.  Surface  and  subsurface  stresses  resulting  from  rubbing  friction  contact  combined 
with  the  stresses  resulting  from  the  overall  fluctuating  stress  field. 

There  are  many  practical  examples  of  structures  in  which  fretting  occurs.  These  structures 
include:  bolted  and  riveted  joints,  suspension  cables,  metallic  coil  ropes,  key-way  shaft  cou¬ 
plings,  coil  wedges  of  generator  rotors  and  blade  dovetail  contact  sections  of  turbine  engines 
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[2].  The  fretting  problem  is  made  worse  by  the  fact  that  many  of  the  structures  that  undergo 
fretting  cannot  be  easily  inspected.  The  first  detection  of  a  crack  may  be  when  the  structure 
fails. 

The  first  reported  observation  of  the  occurrence  of  fretting  during  fatigue  testing  was  made 
by  Eden,  Rose,  and  Cunningham  [3]  in  1911.  They  observed  corrosion  in  the  holders  of  the 
fatigue  test  specimens.  The  corrosion  was  such  that  they  had  difficulty  removing  the  broken 
test  specimens  from  the  testing  apparatus  for  the  heavily  loaded  tests  which  lasted  for  several 
days.  Eden,  Rose,  and  Cunningham  [3]  also  recorded  the  first  attempt  of  the  use  of  palliatives 
for  fretting  fatigue.  They  oiled  the  fatigue  specimens  before  putting  them  in  place,  however 
they  reported  that  this  was  only  partially  effective.  The  specimens  wpre  easier  to  remove,  but 
corrosion  products  were  still  produced.  The  tests  conducted  by  Eden  et  al  [3]  were  not  fretting 
fatigue  tests,  the  observed  fretting  was  noted  but  not  analyzed. 

According  to  Hoeppner  [4]  the  first  systematic  investigation  of  fretting  was  performed  by 
Tomlinson  in  1927,  in  which  he  determined  that  surface  damage  was  related  to  the  relative 
displacement  between  the  contacting  surfaces.  Warlow-Davies  [5]  investigated  the  impact  of 
fretting  on  the  fatigue  properties  of  a  material  in  1941  and  demonstrated  that  fretting  reduced 
fatigue  strength.  He  did  not  perform  fretting  and  fatigue  tests  simultaneously,  but  fretted  the 
specimens  and  then  performed  fatigue  tests  on  the  fretted  specimens.  Warlow-Davies  performed 
the  tests  on  medium-carbon  steel  and  for  a  nickel-cliromium-molybdenum  alloy  steel.  The  loss 
in  fatigue  strength  after  fretting  was  found  to  be  13  percent  and  18  percent  respectively.  By 
separating  the  fretting  and  fatigue  portions  of  the  test,  Warlow-Davies  was  able  to  measure 
the  reduction  in  fatigue  strength  due  only  to  the  surface  damage  caused  by  fretting.  Later 
research,  such  as  that  conducted  by  Waterhouse  [29],  showed  that  simultaneous  fretting  and 
fatigue  resulted  in  much  higher  reductions  in  fatigue  strength.  There  have  been  numerous 
other  experiments  demonstrating  the  deleterious  effects  of  fretting  on  fatigue  strength.  Birch 
[11],  Szolwinski  fc  Farris  [30]  and  Aclibnazari  &  Hoeppner  [15]  are  just  a  few  of  the  other  works 
that  have  also  demonstrated  the  damaging  effects  of  fretting  on  fatigue  life. 

Much  work  has  been  performed  in  fretting  fatigue  since  the  investigation  conducted  by 
Warlow-Davies.  This  thesis  will  not  attempt  to  reference  all  the  major  work  performed  in 
the  area  of  fretting  fatigue.  Hoeppner  in  [4]  provides  a  historical  review  on  fretting  fatigue 
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Figure  1-1:  Simplified  schematic  of  bridge  type  fretting  pads. 

research  as  well  as  a  brief  review  of  the  current  state  of  knowledge,  which  is  part  of  an  entire 
work  devoted  to  the  topic  of  fretting  fatigue.  Birch  [11]  also  provides  a  brief  summary  of 
fretting  research. 


1.2  Motivation  and  Objectives 

Although  the  damaging  effects  of  fretting  on  the  fatigue  properties  of  a  material  have  been 
reported  quite  extensively  since  the  work  of  Warlow- Davies,  there  is  still  not  a  consensus  on 
the  critical  parameter(s)  for  fretting  fatigue.  Waterhouse  [9]  states  that  the  main  factor  in 
producing  fatigue  failure  is  the  additional  shear  stress  in  the  contact  region  caused  by  the 
frictional  forces  between  the  contacting  bodies,  while  Nishioka  and  Hirakawa  [10]  stress  the 
importance  of  the  relative  slip  amplitude  and  Adibnazari  and  Hoeppner  [15]  demonstrate  the 
importance  of  the  normal  pressure  on  fretting  fatigue  life.  According  to  Birch  [11],  Beard  lists 
more  than  fifty  parameters  of  concern  for  fretting  fatigue.  The  lack  of  consensus  on  the  critical 
fretting  parameter(s)  is  likely  due  to  the  experimental  configurations  that  have  been  used  in 
the  past.  Many  of  the  past  fretting  fatigue  experiments  have  been  performed  using  bridge  type 
pads  as  shown  in  Figure  1-1. 

The  bridge  pads  produce  stress  singularities  at  the  contact  edges  making  a  stress  analysis 
very  difficult.  Many  of  the  experimental  configurations  in  the  past  also  lacked  the  ability 
to  measure  all  of  the  fretting  parameters,  specifically  the  oscillating  tangential  contact  load, 
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which  is  a  very  important  fretting  parameter.  If  the  oscillating  tangential  load  is  not  present, 
the  problem  becomes  an  indentation  problem  and  not  fretting.  This  inability  to  measure  the 
contact  load  conditions  again  limited  the  information  available  for  analysis.  The  first  objective 
of  this  work  is  to  determine  the  critical  fretting  fatigue  parameter(s) . 

The  limitations  discussed  above  are  not  in  this  work.  The  work  of  Birch  [11]  and  the  ongoing 
experimental  work  of  Conner  [12]  provide  the  experimental  data  used  in  this  work.  The  testing 
apparatus  used  by  Birch  and  Conner  employed  spherical  fretting  pads.  The  spherical  fretting 
pads  allow  analytical  expressions  for  the  stress  field  due  to  contact  to  be  determined,  which 
was  first  done  by  Hamilton  and  Goodman  [13].  The  testing  apparatus  used  by  Birch  and 
Conner,  described  in  detail  in  [25],  allows  measurement  of  all  the  fretting  parameters  required 
to  determine  the  stress  fields  due  to  contact.  The  expressions  given  in  [13]  were  expressed  again 
later  by  Hamilton  [14]  in  a  more  convenient  explicit  form.  The  form  developed  by  Hamilton 
[14]  is  implemented  in  this  work.  Hamilton  and  Goodman  provide  the  elastic  stress  fields. 
This  does  not  present  a  problem  in  this  work  as  the  stresses  remained  in  the  elastic  region  for 
all  but  a  few  tests. 

The  expressions  given  by  Hamilton  [14]  are  for  a  global  sliding  contact,  which  is  not  the 
condition  for  most  of  the  experiments  analyzed  in  this  work.  Most  of  the  tests  analyzed  in 
this  thesis  have  slip/stick  conditions.  Chivers  and  Gordelier  [16]  show  how  the  Hamilton 
expressions  must  be  modified  to  accommodate  slip/stick  fretting  conditions.  This  modification 
has  been  implemented  in  tliis  work  for  the  slip/stick  experimental  conditions  when  necessary. 
It  will  be  demonstrated  in  this  work  that  the  critical  fretting  parameters  can  be  captured  by 
one  parameter,  the  maximum  von  Mises  effective  stress  at  the  surface  due  to  contact. 

The  second  objective  of  this  work  is  to  analyze  fretting  fatigue  life  prediction  methods 
using  past  experimental  work.  A  widely  accepted  fretting  fatigue  life  prediction  method  does 
not  yet  exist.  The  stress  life  approach,  which  is  widely  used  in  industry  today,  is  examined 
for  its  applicability  for  fretted  components.  There  have  been  several  life  prediction  methods 
proposed  for  fretting  fatigue,  but  none  have  utilized  the  stress-life  approach.  Edwards  [6]  in 
1981  proposed  a  fretting  fatigue  life  prediction  method  using  linear  elastic  fracture  mechanics. 
Faanes  and  Fernando  [7]  in  1994  also  proposed  a  life  prediction  method  using  linear  elastic 
fracture  mechanics,  but  added  a  short  crack  correction  to  account  for  the  effect  of  the  high 


stress  gradient  due  to  fretting  on  short  crack  growth.  The  success  of  linear  elastic  fracture 
mechanics  applied  to  fretting  has  been  limited  by  the  difficulty  in  developing  stress  intensity 
factors  to  be  used  to  predict  short  crack  growth  rates.  Giannakopoulos,  Lindley  and  Suresh 
[46]  propose  a  method  which  circumvents  tliis  problem  for  sharp  edged  contacts  in  their  crack 
analogue  model.  Szolwinski  and  Farris  [30]  in  a  recent  paper  propose  a  fretting  fatigue  life 
prediction  method  that  uses  the  Smith- Watson-Topper  expression  [8]  to  predict  the  number  of 
cycles  to  crack  nucleation  and  employs  linear  elastic  fracture  mechanics  to  predict  the  long  crack 
growth  rate,  thereby  avoiding  the  short  crack  growth  rate  problem.  This  work  will  examine 
both  the  stress-life  approach  and  the  damage  tolerant  approach  used  to  predict  fatigue  life. 

The  final  portion  of  this  work  analyzes  fretting  fatigue  palliatives;  coatings  and  shot  peening. 
Shot  peening  has  been  shown  by  Waterhouse  [38],  Bignonnet  [39],  and  Olivers  and  Gordelier 
[41]  to  be  an  effective  fretting  fatigue  palliative.  However,  a  detailed  analysis  of  the  shot  peening 
process  applied  to  fretted  components  with  the  objective  of  determining  the  optimum  peening 
depth  has  not  yet  been  performed.  Fretting  produces  large  tensile  and  compressive  stresses  at 
the  contact  edges  and  both  must  be  considered  when  determining  the  optimum  peening  depth. 
A  systematic  method  to  determine  the  optimum  peening  depth  for  a  fretted  component  is  given 
in  this  work,  using  the  elastic  stress  field  expressions  provided  by  Hamilton.  This  analysis  also 
includes  the  effects  of  changing  the  coefficient  of  friction  on  the  optimal  peening  depth.  The 
elastic  stress  fields  due  to  contact  will  also  be  used  to  provide  a  brief  discussion  of  the  benefits 
of  protective  coatings  for  fretting  fatigue. 

This  thesis  is  organized  in  the  following  manner:  Chapter  2  provides  a  description  of  the 
sphere  on  flat  plane  contact  mechanics  used  in  this  work.  Chapter  3  uses  the  works  of  Hamilton 
[14],  and  Chivers  and  Gordelier  [16]  to  analyze  the  experimental  results  provided  by  Birch  [11], 
and  Conner  [12].  Chapter  4  analyzes  fretting  fatigue  life  prediction.  Chapter  5  discusses 
fretting  fatigue  palliatives  and  provides  a  method  to  determine  the  optimum  peening  depth  for 
a  fretted  component.  Chapter  6  provides  the  conclusions  from  tliis  work  and  suggestions  for 
further  study  in  work  related  to  fretting  fatigue. 
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Chapter  2 


Review  of  Sphere  on  Flat  Plane 
Contact  Mechanics 


2.1  Introduction 

Before  the  analysis  of  Hamilton  [14]  and  Olivers  and  Gordelier  [16]  is  applied  to  experimental 
tests,  a  brief  review  of  the  sphere  on  flat  plane  contact  mechanics  is  presented.  The  expressions 
reviewed  here  serve  as  the  foundation  for  the  analysis  performed  in  this  work.  The  discussion 
begins  with  the  analysis  conducted  by  Hertz  in  1882. 


2.2  Hertzian  Contact 


Over  one  hundred  years  ago  Hertz  analyzed  the  problem  of  normal,  frictionless  contacting 
spheres.  Hertz  solved  for  the  contact  pressure  distribution,  p(r)  and  the  circular  radius  of 
contact,  a,  by  solving  the  integral  equation  governing  the  boundary  conditions.  The  radius  of 
contact  as  determined  by  Hertz  [17]  is  given  by: 


f  SirPRiRi 

\mTW) 


T-t/f  | 
nEi  E*2 


(2.1) 


where  E  and  v  are  Young’ s  Modulus  and  Poisson’s  ratio  respectively  and  R  is  the  radius  of 
curvature.  The  indices  1,2  refer  to  the  two  contacting  bodies.  The  solution  derived  by  Hertz 


Figure  2-1:  (a)  Fretting  contact  between  a  shpere  and  a  fiat  surface,  (b)  Top  view  showing 
circular  contact  region  and  stick  zone,  c. 

for  normal  contacting  spheres  can  be  applied  to  a  sphere  on  a  flat  plane  geometry  by  noting 
that  =  R^~  +  R^-e>  Rpiane  =  oo,  therefore  R  =  Reff  =  R3phere>  Huber  [18]  in  1904, 
was  the  first  to  solve  for  the  stress  fields  beneath  the  contacting  surfaces. 


2.2.1  Sliding  Elastic  Sphere  on  a  Flat  Surface 

The  indentation  problem  solved  by  Hertz  and  Huber  becomes  a  fretting  problem  when  an 
oscillatory  tangential  force  is  applied  to  the  sphere.  The  problem  geometry  is  shown  in  Figure 
2-1. 

Once  the  tangential  force,  Q,  is  introduced,  sliding  or  partial  slip  is  also  introduced.  The 
Amonton-Coulomb  law  relates  the  normal  pressure  and  the  tangential  tractions,  q  by: 


M  _  M  _  „ 

P  ~  p 


(2.2) 


where  ji  is  the  coefficient  of  kinetic  friction. 

As  given  first  by  Mindlin  [19]  in  1949  and  again  later  by  Hamilton  [14]  the  stresses  at  the 


contact  surface  for  incipient  sliding  conditions  are  given  by: 


o-r  —  Va?  -  r2,  0  <  r  <  a  (2.3) 

2tt  cr 

\]al  -r'1,  0  <  r  <  a  (2.4) 

ry2  and  the  boundary  stresses  outside  the  circular  area  of  contact  are  zero.  As  mentioned 
above,  these  boundary  conditions  only  apply  under  conditions  of  gross  sliding,  that  is  when 
Q  =  fiP.  The  expression  given  by  Equation  2.3  is  also  the  expression  for  contact  pressure  due 
to  indentation  without  sliding  and  the  expression  given  for  the  contact  radius  by  Equation  2.1 
is  also  unchanged  by  the  addition  of  a  tangential  load.  Therefore,  the  presence  of  the  shear 
stresses  on  the  surface  has  not  affected  either  the  circular  contact  radius  or  the  contact  pressure. 


2.2.2  Partial  Sliding 

When  the  applied  tangential  force  is  less  than  the  product  of  fiP \  macroscopic  sliding  does  not 
occur,  however  frictional  traction  at  the  contact  interface  leads  to  small  relative  displacement 
over  part  of  the  interface.  This  small  relative  motion  is  commonly  referred  to  as  micro-slip. 
The  remaining  portion  of  the  interface  deforms  without  any  relative  motion.  The  region  in 
which  no  relative  motion  occurs  is  referred  to  as  the  stick  region.  All  the  points  on  the  surface 
within  the  stick  region  undergo  the  same  tangential  displacements.  The  stick  region  is  shown 
in  Figure  2-1.  Mindlin  [20]  has  shown  that  a  stick  region  exists  in  the  contact  interface  and 
the  radius,  e,  of  the  outer  boundary  of  the  stick  region  to  be  given  by: 


c 

a 


1  - 


?p\ 


1 1/3 


c  <r  <  a. 


(2.5) 


As  shown  in  Figure  2-1,  the  stick  region  is  concentric  with  the  contact  area,  this  geometric 
relationship  is  also  observed  for  experimental  specimens. 

Just  as  the  boundary  conditions  had  to  be  reestablished  for  sliding,  they  must  again  be 
modified  to  account  for  micro-slip.  The  stresses  at  the  contact  surface  are  now  given  by  [19]: 


7~zx  — 


3  jiP 

27 ra2 


\J a 2  —  r2,  c  <r  <  a  and 


(2.6) 
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(2.7) 


The  expression  for  the  normal  pressure  is  the  same  as  in  the  global  sliding  condition.  The 
relative  tangential  displacements  are  given  by  [19]: 


r  _  3//P  f  2  -  vi  2  -  Vj 
~  16a  \  Gi  +  G2 


2/3 > 


(2.8) 


2.3  Hamilton’s  Analysis 

As  mentioned  in  the  previous  section,  the  stress  fields  beneath  the  surface  of  a  circular  indenter 
were  first  derived  by  Huber  in  1904.  However,  the  elastic  stress  fields  beneath  the  surface  of  a 
circular  sliding  contact  on  a  flat  surface  were  not  derived  until  1966  by  Hamilton  and  Goodman 
[13].  Mindlin’s  work  in  circular  sliding  contact  focused  primarily  on  surface  compliance  and 
the  micro-slip  problem.  He  did  not  analyze  the  stresses  below  the  contact  surface. 

It  must  be  stressed  that  the  equations  given  by  Hamilton  are  for  a  global  sliding  sphere 
and  not  for  stick /slip  conditions.  It  was  shown  earlier  that  the  boundary  stresses  are  not 
the  same  for  the  two  conditions.  The  stress  equations  for  the  global  sliding  condition  must 
be  modified  slightly  to  obtain  the  correct  stresses  for  situations  involving  micro-slip.  Chivers 
and  Gordelier  [16]  demonstrate  how  the  expressions  must  be  modified.  The  results  of  a  stress 
analysis  performed  using  the  expressions  given  by  Hamilton  are  shown  in  Figure  2-2.  The 
expressions  used  to  generate  this  plot  are  given  in  Appendix  D.  The  material  used  to  generate 
this  plot  was  A1  7075-T6.  Some  material  properties  for  this  alloy  are  listed  in  Table  2.1  [11] 
[21]  [22].  Spherical  fretting  pads  with  a  one  inch  radius  were  used  and  global  sliding  is  assumed. 
Using  the  values  for  fi)  P,  and  v  in  Table  2.1  and  the  elastic  stress  field  equations  given  by 
Hamilton,  the  stresses  in  the  x  and  y  directions  have  been  plotted  in  Figure  2-2  along  with 
the  von  Mises  effective  stress,  creff ,  at  2  =  0,  y  =  0.  The  stresses  have  been  normalized  by 
the  maximum  contact  pressure,  Po>  Po  =  The  von  Mises  effective  stress  parameter,  crcff, 

has  been  used  because  this  is  a  multiaxial  problem.  The  effective  stress  is  used  to  convert  the 
multiaxial  stress  state  to  a  single  stress  parameter.  The  expression  for  the  effective  stress  is 

IS 


Material  Property 

Value 

Young’s  modulus  (E) 

72  GPa 

Poisson  ratio  ( v ) 

0.3 

Monotonic  yield  stress  (cry) 

469  MPa 

Coefficient  of  kinetic  friction  (p) 

1.2 

Tensile  strength  (ots) 

579  MPa 

Endurance  limit  (crend) 

179  MPa 

Fatigue  strength  coefficient  (o') 

1317  MPa 

Fatigue  strength  exponent  ( b ) 

-0.126 

Fatigue  limit  stress  (<7d) 

230  MPa 

Fatigue  ductility  coefficient  (ef) 

0.19 

Fatigue  ductility  exponent  (c) 

-0.56 

Table  2.1:  Properties  of  7075-T6  aluminum  alloy 


given  by: 


<Tcft  — 


(o-x  -  <7y)2  +  (cry  -  crz)2  +  (er2  -  crx)2  +  6  (r2y  +  T*x  +  r*x) 


1/2 


(2.9) 


The  stress  fields  displayed  in  Figure  2-2  are  not  the  only  non-zero  stress  fields  at  the  surface; 
the  shear  stress  rzx  and  the  normal  stresses  in  the  z-direction  are  also  non-zero,  however  they 
both  go  to  zero  at  the  contact  edges.  The  stresses  below  the  surface  at  a  depth  of  20  fim  is 
plotted  in  Figure  2-3.  Comparing  Figures  2-2  and  2-3,  it  can  be  seen  that  the  stresses  due  to 
contact  quickly  decay  away,  especially  the  tensile  stress  at  the  trailing  edge  of  contact,  x  =  —a, 
in  the  x-direction.  From  Figure  2-2  it  can  be  seen  that  there  is  a  large  peak  in  crx  at  the  trailing 
edge  of  contact.  This  stress  is  the  highest  stress  on  the  plot.  In  Figure  2-3,  the  peak  tensile 
stress  at  the  trailing  edge  is  much  lower  in  magnitude  and  the  slope  of  the  <7X  curve  is  also  much 
smoother. 

The  maximum  von  Mises  effective  stress  is  also  at  the  surface.  This  is  true  for  a  sphere 
sliding  on  a  flat  surface  when  (i  >  0.3.  When  fi  =  0,  or  when  the  tangential  force,  Q,  is  zero, 
the  maximum  von  Mises  stress  shifts  below  the  surface  to  the  point  z  =  0.5a,  for  v  =  0.3  [14]. 
However,  for  a  sphere  fretting  on  a  flat  plane,  the  maximum  von  Mises  and  the  maximum  tensile 
stresses  generally  occur  at  the  surface.  Although  the  maximum  tensile  stress  has  a  severe  peak 
at  the  trailing  edge  of  contact  at  the  surface,  the  von  Mises  effective  stress  is  roughly  constant 
over  a  very  large  region  as  can  be  seen  in  Figure  2-2. 


IQ 


The  analysis  performed  in  this  work  will  focus  on  the  surface  stresses,  at  the  trailing  edge  of 
contact,  with  the  exception  of  the  shot  peening  analysis,  which  analyzes  below  the  surface  and 
inside  the  contact  edges.  The  contact  stresses  are  maximum  at  the  surface  and  the  maximum 
tensile  stress  occurs  in  the  material  at  the  trailing  edge  of  the  contact;  therefore,  this  work 
will  focus  its  analysis  on  this  location.  The  von  Mises  effective  stress  is  just  as  high  at  other 
locations  on  the  surface,  however  this  is  in  general  due  to  large  compressive  stresses  at  that 
location.  It  is  assumed  that  the  large  tensile  stress  is  responsible  for  crack  nucleation.  The 
contact  stresses  given  by  Hamilton  can  be  greatly  simplified  for  x  =  -a,  z  =  0  and  y  —  0,  the 
problem  becomes  a  biaxial  one.  The  only  non-zero  stresses  at  this  location  are  ctx  and  cry,  and 


are  given  by: 


a*  =  Po 


(1  —  2j/)  4  +  t'  , 

i — - — -  ±  — 3 — ,  x  =  Ta,  and 

O  o 


<Jy  =  Po 


(2  v  —  1)  Z-Kvy, 

3  8 


CO  — *  -pfl. 


(2.10) 


(2.11) 


Note  carefully  the  position  of  the  signs.  The  maximum  tensile  stress  in  both  the  x  and  y 
directions  occurs  at  the  trailing  edge  of  contact.  Again,  it  must  be  emphasized  that  the 
stresses  given  in  equations  2.10  and  2.11  are  for  full  sliding  and  not  for  stick/slip  conditions. 
Figure  2-4  displays  the  curves  for  ax  under  incipient  sliding  and  stick/slip  conditions.  As  can 
be  seen  in  this  figure,  when  stick/slip  is  accounted  for,  crx  is  reduced  for  the  given  contact  load 
conditions.  This  figure  was  generated  using  a  c/a  ratio  of  0.5. 

Recall  that  the  ratio  of  c/a  indicates  how  close  to  full  sliding  or  full  stick  the  fretting 
condition  is.  As  c/a  0,  global  sliding  is  introduced,  when  c/a  =  1,  full  stick  conditions 
exist.  From  Figure  2-4,  it  can  be  seen  that  the  magnitude  of  the  maximum  tensile  stress 
is  reduced  when  micro-slip  is  taken  into  consideration.  This  fact  will  become  more  relevant 
later  during  the  shot-peening  analysis  of  this  work.  The  reduction  in  stress  under  stick/slip 
conditions  is  in  direct  conflict  with  experimental  results.  Experiments  have  shown  that  fretting 
is  most  damaging  to  the  fatigue  properties  of  a  system  under  micro-slip  conditions  with  relative 
displacements  of  10-50  /im.  The  introduction  of  global  sliding  improves  the  fretting  fatigue 
life.  It  has  been  postulated  that  this  improvement  under  sliding  conditions  may  be  due  to 
the  wearing  away  of  fatigue  cracks  from  the  surface,  so  that  they  do  not  become  propagating 
fatigue  cracks.  Tlxis  is  one  example  of  the  limitations  of  the  stress  analysis. 
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x/a 

Figure  2-4:  Plot  demonstrating  the  effects  of  stick-slip  conditions  on  crx  for  A1  7075-T6  (c/a  = 
0.5). 

In  most  practical  structures  and  scientific  experiments,  the  tangential  load,  Q,  which  changes 
the  problem  from  one  of  indentation  to  fretting,  is  caused  by  an  applied  cyclic  bulk  stress. 
Therefore,  to  perform  a  complete  analysis,  the  stresses  induced  in  the  material  by  the  bulk  load 
must  also  be  accounted  for.  If  the  material  behaves  elastically,  the  principle  of  superposition 
can  be  applied.  Using  superposition,  the  stresses  at  the  surface,  at  r  =  a,  are  given  by: 

ffxtot  =  <?xc  +  o"xb  and  (2.12) 

U"  y  tOt  =  0"yC  +  Cyb,  (2.13) 

where  crxtot  and  crytol  are  the  total  stresses  in  their  respective  directions.  crxc  and  ayc  are 
the  contact  stresses  in  their  respective  directions  and  crxb,  0yb  are  the  bulk  stresses  in  the  x 
and  y  directions  respectively.  The  condition  described  above,  that  of  a  cyclic  bulk  load  in 
combination  with  the  applied  contact  loads,  was  the  method  used  to  produce  the  experimental 
results  analyzed  in  this  work.  The  bulk  load  also  produces  eccentricity  in  the  contact  geometry. 
The  contact  area  and  stick  zone  are  no  longer  perfectly  circular.  This  effect  is  discussed  in 
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detail  in  [23]. 
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Chapter  3 


Experimental  Results 


3.1  Introduction 

As  shown  earlier,  the  elastic  stress  fields  for  a  sphere  sliding  on  a  flat  plane  are  known  and  can 
be  found  as  functions  of  the  normal  and  tangential  loads.  The  testing  apparatus  developed  at 
the  Massachusetts  Institute  of  Technology  gives  the  user  the  ability  to  control  and  measure  both 
the  normal  and  tangential  contact  loads.  This  apparatus  has  been  used  to  conduct  fretting 
fatigue  tests  on  A1  7075-T6  and  Ti-6A1-4V.  The  testing  apparatus  is  described  in  [25].  A 
simplified  diagram  of  the  test  setup  is  shown  in  Figure  3-1.  These  tests  provide  the  necessary 
information  to  calculate  the  elastic  stress  fields  due  to  contact,  which  will  be  discussed  in  detail 
in  the  following  sections.  The  elastic  stress  fields  at  the  surface  due  to  contact  do  not  completely 
describe  the  fretting  phenomenon,  there  are  some  shortcomings.  The  stress  analysis  does  not 
account  for  the  condition  of  the  surface  and  may  not  give  the  correct  location  of  crack  initiation. 
However,  the  stress  analysis  provides  a  very  good  starting  point  for  a  systematic  analysis  of  the 
problem. 

3.2  A1  7075- T6  Experimental  Results 

The  experimental  data  provided  in  the  work  performed  by  Birch  [11]  and  Conner  [12]  is  shown 
in  Table  3.1.  Table  3.1  provides  all  the  necessary'  information  to  calculate  the  stress  fields.  The 
bulk  stress,  <r&,  for  these  experiments  was  applied  in  the  x-direction  only,  with  a  stress  ratio, 


9  A 


Applied  Bulk  Load 


Figure  3-1:  Simplified  diagram  of  the  fretting  fatigue  testing  apparatus. 


R=- 1. 

Care  must  be  taken  when  superposing  the  bulk  stress.  As  discussed  earlier  the  trailing 
edge  of  contact  is  the  point  of  interest  for  this  analysis;  however,  during  each  complete  stress 
cycle  the  trailing  edge  becomes  the  leading  edge  as  the  bulk  stress  is  reversed.  The  trailing 
edge  during  the  tensile  portion  of  the  bulk  stress  loading  cycle  becomes  the  leading  edge  during 
the  compressive  portion  of  the  loading  cycle.  This  location  on  the  specimen  undergoes  the 
highest  tensile  stresses  and  the  largest  stress  amplitude.  The  opposite  edge  of  contact  endures 
lower  stresses  because  it  becomes  the  trailing  edge  of  contact  during  the  compressive  portion 
of  the  bulk  loading  cycle.  It  would  be  expected  that  the  cracks  would  initiate  from  the  edge 
of  contact  wrhich  undergoes  the  highest  stresses.  This  fact  is  verified  in  [23]  which  displays  a 
fretting  scar  with  cracking  visible  only  on  one  side  of  the  circular  contact  region. 

The  loading  conditions  displayed  in  Table  3.1  were  used  to  calculate  the  surface  stresses 
at  the  trailing  edge  of  contact  for  each  test  specimen.  To  better  analyze  the  results  of  the 
stress  analysis,  this  work  has  isolated  two  series  of  experiments,  which  correlate  to  two  different 
loading  conditions.  The  first  series  of  tests  includes  those  tests  where  the  bulk  stress  was 
roughly  constant  (^83MPa),  while  the  contact  loads  were  varied.  Therefore,  the  effects  of  the 


Test 


Dest 

P  (N) 

Q(N) 

Q/(m  x-P) 

o-b  (MPa) 

Nf  (Cycles) 

20.8 

15 

0.60 

83 

449,5500 

30 

15 

0.42 

85 

480,000 

15.6 

15 

0.80 

85 

395,000 

12.5 

15 

1.0 

83 

361,000 

5 

18.5 

13.5 

0.61 

77 

551,000 

6 

16 

11.7 

0.61 

83 

530,000 

7 

8.75 

7 

0.67 

83 

3,300,000  (no  failure) 

8 

13 

7 

0.45 

83 

10,680,000  (no  failure) 

9 

13.9 

10 

0.60 

83 

803,000 

10 

10.3 

7.5 

0..61 

83 

2,940,000 

11 

7.31 

6.6 

0.75 

83 

3,450,000  (no  failure) 

12 

8.33 

8 

0.80 

83 

616,000 

13 

20 

15 

0.63 

69 

1,410,000  (no  failure) 

14 

20 

21.5 

0.89 

54 

2,186,000  (no  failure) 

15 

20 

20 

0.83 

61 

2,500,000  (no  failure) 

16 

24 

20.25 

0.70 

61 

2,666,000  (no  failure) 

17 

28 

20 

0.59 

61 

2,505,000  (no  failure) 

18 

11.5 

8.6 

0.62 

83 

2,608,000  (no  failure) 

19 

20 

15 

0.63 

83 

549,000 

20 

20 

15 

0.63 

70 

516,000 

21 

20 

15 

0.63 

56 

1,540,000  (no  failure) 

22 

20 

15 

0.63 

63 

2,940,000  (no  failure) 

23 

20 

15 

0.63 

59 

1,777,000  (no  failure) 

24 

25 

23 

0.77 

270 

297,000 

25 

25 

22.5 

0.75 

270 

135,000 

26 

22.5 

19.5 

0.72 

230 

480,000 

27 

20 

17 

0.71 

195 

575,000 

28 

20 

15 

0.63 

75 

3,050,000 

29 

20 

15 

0.63 

100 

390,000 

30 

20 

15 

0.63 

125 

258,063 

31 

20 

15 

0.63 

150 

181,118 

32 

20 

15 

0.63 

74 

2,497,226 

33 

20 

15 

0.63 

67 

3,200,845  (no  failure) 

Table  3.1:  A1  7075-T6  fretting  fatigue  experiments 


Test 

a  (pm) 

c/a 

o’  xtot  (MPa) 

<Tec  (MPa) 

Nf  (Cycles) 

8 

184 

0.82 

338" 

248 

10,680,000  (no  failure) 

7 

162 

0.69 

360 

264 

3,300,000  (no  failure) 

11 

152 

0.63 

362 

264 

3,450,000  (no  failure) 

10 

171 

0.73 

361 

266 

2,940,000 

18 

177 

0.72 

375 

2,608,000  (no  failure) 

12 

159 

0.58 

384 

285 

616,500 

9 

188 

0.74 

388 

292 

803,000 

6 

197 

0.73 

530,000 

2 

244 

0.83 

317 

480,000 

1 

216 

0.74 

432 

334 

449,500 

19 

213 

1 

337 

549,000 

3 

196 

s 

456 

351 

395,000 

4 

182 

473 

366 

361,000 

Table  3.2:  A1  7075-T6  Constant  Bulk  load  data  analysis 


bulk  stress  are  the  same  in  each  test  and  the  difference  in  life  is  due  only  to  the  difference  in 
contact  stress.  The  second  series  of  tests  analyzed  are  made  up  of  tests  with  constant  contact 
loading  (P  =  20  N,  Q  =  15  N),  wdiile  the  bulk  stress  is  varied,  which  isolated  the  effect  of  the 
bulk  stress  on  fatigue  life. 

The  results  of  the  stress  calculations  for  the  constant  ai  series  of  tests  are  shown  in  Table 
3.2.  Table  3.2  displays  the  calculated  values  for  the  contact  radius,  a,  the  ratio  of  the  stick 
zone  radius  to  the  contact  zone  radius  (c/a),  the  maximum  tensile  stress  at  the  trailing  edge  of 
contact,  <7 xtot  and  the  effective  stress  due  to  contact  at  the  trailing  edge,  <jec. 

Table  3.2  has  been  arranged  in  order  from  lowest  to  highest  <rec.  Note  that  the  fatigue 
lifetimes  as  a  result  are  nearly  arranged  in  order  from  greatest  number  of  cycles  to  lowest. 
This  fact  indicates  that  a  relationship  exists  between  the  surface  stresses  at  the  trailing  edge  of 
contact  and  the  fretting  fatigue  life  of  the  specimen.  This  relationship  will  be  discussed  further 
in  the  next  chapter,  where  fretting  fatigue  life  prediction  methods  are  discussed.  The  fatigue 
life  versus  effective  stress  due  to  contact  is  plotted  in  Figure  3-2.  Figure  3-2  again  demonstrates 
the  existence  of  a  relationship  between  the  fatigue  life  and  the  surface  stresses  due  to  contact. 
Recall  that  each  one  of  these  tests  were  conducted  with  approximately  the  same  applied  bulk 
load,  therefore  the  difference  in  fatigue  life  is  due  to  the  difference  in  contact  loading  only.  It 
can  be  seen  that  as  the  contact  stresses  are  reduced,  the  fatigue  life  increases.  For  this  bulk 
loading  condition  it  appears  that  specimen  failure  will  not  occur  if  aec  is  maintained  less  than 
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Figure  3-2:  Plot  of  effective  stress  due  to  contact  versus  number  of  cycles  to  failure  for  crb  ~  83 
MPa. 

about  264  MPa,  which  provides  an  endurance  type  stress  combination.  Also  note  that  in 
none  of  the  tests  was  the  effective  stress  due  to  contact  greater  than  the  yield  strength  of  the 
material;  however,  in  test  four,  the  total  stress  in  the  x-direction  does  slightly  exceed  the  yield 
strength  of  the  material. 

The  stress  values  are  not  the  only  important  values  displayed  in  Table  3.2;  the  values 
calculated  for  a  and  c/a  are  also  quite  important.  They  provide  analytical  values  that  can  be 
easily  checked  against  the  actual  scars  on  the  specimens,  which  is  done  later  in  this  chapter. 
This  gives  insight  into  the  validity  of  the  tests. 

The  value  of  c/a  not  only  provides  a  means  to  verify  testing  procedures,  it  also  is  a  key 
parameter  in  characterizing  the  nature  of  the  fretting  condition.  Recall  that  as  c/a  goes  to 
zero,  global  sliding  is  introduced  and  when  c/a  ==  1,  complete  stick  is  occurring.  It  can  be 
seen  from  Table  3.2,  that  in  Test  4  global  sliding  should  be  occurring.  The  c/a  value  can  also 
provide  insight  into  the  existence  of  a  normal  pressure  threshold  discussed  by  Adibnazari  and 
Hoeppner  [15].  Table  3.3  displays  tests  with  the  same  value  of  tangential  load  (15N),  with 
increasing  values  of  normal  load.  It  can  be  seen  that  as  the  normal  load  is  increased,  the 
maximum  tensile  stress  and  the  effective  stress  due  to  contact  actually  decrease.  Therefore,  for 


Test 

P  (N) 

Q(N) 

0-xtot  (MPa) 

<rec  (MPa) 

4 

12.5 

mm 

473 

366 

3 

15.6 

mm 

456 

351 

14 

20 

15 

435 

337 

1 

20.8 

15 

432 

334 

2 

30 

15 

409 

317 

Table  3.3:  Stress  variation  as  the  normal  load  changes  for  a  constant  maximum  tangential  load. 

a  given  tangential  load,  as  the  normal  load  is  increased,  the  contact  stresses  at  the  trailing  edge 
of  contact  decrease,  and  from  Table  3.2  and  Figure  3-2  the  fretting  fatigue  life  should  increase. 
This  provides  a  rational  explanation  of  why  a  pressure  threshold  exists.  The  stresses  decrease 
as  the  normal  load  is  increased  due  to  the  lowering  of  c/a.  As  c/a  goes  to  zero,  the  relative 
displacements  also  go  to  zero  and  complete  stick  results. 

The  results  for  the  second  series  of  tests,  in  which  the  contact  loads  were  maintained  constant 
while  the  bulk  load  was  varied,  are  shown  in  Table  3.4.  The  data  in  Table  3.4  is  presented 
in  order  of  ascending  bulk  load.  The  values  for  contact  radius,  a  and  c/a  are  unchanged  for 
each  test,  since  the  contact  loads  are  identical  for  each  test.  The  values  of  a  and  c/a  were 
found  to  be  213  fin i  and  0.72  respectively.  The  value  of  cja  was  <1,  therefore  all  the  tests 
were  conducted  under  stick/slip  conditions.  The  most  interesting  aspect  of  this  set  of  tests 
is  the  appearance  of  a  fretting  fatigue  endurance  type  limit  for  the  given  contact  conditions. 
It  appears  that  when  the  bulk  stress  is  reduced  below  approximately  70  MPa  for  the  given 
contact  loading  conditions,  specimen  failure  will  not  result.  Figure  3-3  is  a  plot  of  the  bulk 
stress  versus  the  number  of  failures  to  cycle  (S-N  curve)  for  these  contact  conditions.  This 
is  not  a  curve  with  universal  applications.  It  only  applies  to  the  contact  conditions  used  in 
this  set  of  tests.  A  different  set  of  contact  load  conditions  will  result  in  a  different  endurance 
strength.  However,  combining  the  data  from  the  first  set  of  A1  7075-T6  tests  with  this  set 
of  data  and  the  known  endurance  strength  of  the  material  a  curve  that  does  have  universal 
applications  can  be  proposed.  This  curve  is  shown  in  Figure  3-4.  In  this  figure  the  so  called 
knock-down  factor,  which  is  the  factor  reduction  of  the  endurance  strength  of  the  material  due 
to  fretting,  is  plotted  versus  the  effective  stress  due  to  contact  at  the  trailing  edge. 

The  curve  displayed  in  Figure  3-4  is  a  very  rough  approximation  of  the  reduction  in  the 
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Test 

a  (/im) 

c/a 

Cxtot  (MPa) 

<rec  (MPa) 

Nf  (Cycles) 

ffb(MPa) 

21 

213 

Egg] 

408 

337 

1,540,000  (no  failure) 

56 

23 

213 

HI 

411 

337 

1,777,000  (no  failure) 

59 

22 

213 

1 m 

415 

337 

2,940,000  (no  failure) 

63 

33 

213 

0.72 

419 

337 

3,200,845  (no  failure) 

67 

13 

213 

0.72 

421 

337 

1,410,000  (no  failure) 

69 

20 

213 

0.72 

422 

337 

516,000 

70 

32 

213 

0.72 

426 

337 

2,497,226 

74 

28 

213 

0.72 

427 

337 

3,050,000 

75 

19 

213 

0.72 

435 

337 

549,000 

83 

29 

213 

0.72 

452 

337 

390,000 

100 

30 

213 

0.72 

477 

337 

258,063 

125 

31 

213 

0.72 

502 

337 

181,118 

150 

Table  3.4:  A1  7075-T6  stress  analysis  for  constant  contact  loading 


Figure  3-3:  Number  of  cycles  to  failure  versus  <7*>,  S  —  N  curve  for  constant  contact  load 
conditions. 


Table  3.5:  Data  used  to  plot  the  knock-down  factor  curve. 

endurance  strength  of  the  material  due  to  fretting.  The  knock-down  factor  is  given  by: 

Kr  _  Cend( no  fretting) 

”  &  end  (fretting) 

where  aend  is  the  endurance  strength  of  the  material.  Figure  3-4  was  constructed  using  only 
three  data  points;  several  more  series  of  tests  at  different  bulk  loads  need  to  be  performed  to 
determine  the  accuracy  of  this  plot.  However,  the  accuracy  of  this  plot  is  not  the  critical  issue, 
the  premise  on  which  the  plot  was  constructed  is  much  more  important.  The  plot  indicates 
that  K F  is  not  a  single  number  for  a  given  material,  but  it  can  be  determined  if  the  contact 
stresses  are  known.  The  fact  that  the  strength  reduction  changes  for  various  contact  loads  can 
easily  be  seen  in  Figure  3-2.  This  idea  will  be  examined  closer  in  the  next  chapter.  The  data 
points  used  to  make  Figure  3-4  are  displayed  in  Table  3.5.  When  the  contact  stresses  go  to 
zero,  the  knock-down  factor  goes  to  one. 

3.3  Ti-6A1-4V  Experiments 

The  second  material  analyzed  in  this  work  is  Ti-6A1-4V  [12] .  The  same  test  apparatus  discussed 
earlier  w'as  also  used  in  the  Ti-6A1-4V  fretting  fatigue  tests,  utilizing  the  same  sphere  on  flat 
surface  fretting  geometry.  Table  3.6  provides  the  necessary  material  properties  to  perform  the 
stress  calculations  for  this  material.  Fuchs  [22]  provides  fatigue  properties  for  many  different 
materials  and  was  used  along  with  [11]  and  [45]  to  produce  Table  3.6, 

The  loading  conditions  used  to  conduct  the  tests  on  the  titanium  alloy  are  presented  in 
Table  3.7.  The  results  of  the  stress  analysis  are  displayed  in  Table  3.8.  The  value  used  for  the 
coefficient  of  friction  was  (i  =  0.5  for  the  displayed  results.  The  fretting  pad  radius,  R,  used  in 
these  tests  was  0.5  inches  versus  one  inch  used  for  the  A1  7075-T6  tests. 

A  plot  of  the  fatigue  life  versus  the  effective  stress  due  to  contact  is  shown  in  Figure  3-5. 
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Table  3.6:  Ti-6A1-4V  material  properties. 


Test 

P  (N) 

Q  (N) 

Q/(M  x  P) 

<rb  (MPa) 

Nf  (Cycles) 

1 

n 

31.5 

225 

691,612 

2 

H 

23.5 

0.94 

400 

109,574 

3 

50 

14.5 

0.58 

394 

107,879 

4 

50 

13 

0.52 

357 

31,540  (interrupted) 

5 

50 

12.25 

0.48 

298 

304,456 

6 

50 

14.5 

0.58 

300 

474,656 

7 

50 

23 

0.96 

327 

143,683 

8 

30 

18 

1.20 

300 

258,576 

9 

30 

9 

0.60 

300 

3,257,000  (no  failure) 

Table  3.7:  Test  parameters  for  Ti-6A1-4V  fretting  fatigue  tests 
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Test 

a  (^m) 

c/a 

<Txtot  (MPa) 

o-ec  (MPa) 

1 

195 

0 

973 

714 

2 

195 

0.39 

988 

567 

3 

195 

0.75 

839 

448 

4 

195 

0.78 

775 

426 

5 

195 

0.80 

702 

414 

6 

195 

0.75 

745 

448 

7 

195 

0.43 

907 

564 

8 

164 

0 

904 

577 

9 

164 

0.74 

682 

384 

Table  3.8:  Ti-6A1-4V  stress  analysis 


This  figure  includes  tests  5-9,  where  the  bulk  stress  was  ~  300  MPa.  Again,  the  fatigue  life 
increases  as  the  contact  stresses  are  reduced.  An  endurance  type  limit  is  also  seen  for  this 
material.  For  at  ~  300  MPa,  it  appears  specimen  failure  will  not  occur  if  <rec  is  maintained 
less  than  about  390  MPa.  More  experimental  tests  are  needed  to  confirm  this  observation. 
Not  enough  experimental  data  was  available  to  plot  a  graph  of  the  fretting  fatigue  life  versus 
applied  bulk  stress. 

3.4  Experimental  Observations 

Observations  of  the  fretting  pads  and  specimens  have  been  made  using  a  scanning  laser  micro¬ 
scope  and  an  optical  microscope.  The  pictures  provided  in  this  section  were  produced  using 
images  provided  by  the  scanning  laser  microscope.  The  scanning  laser  microscope  made  it 
very  simple  and  quick  to  perform  accurate  measurements  of  the  fretting  parameters,  such  as 
the  contact  radius,  a  and  the  radius  of  the  stick  zone,  c. 

The  expressions  for  the  contact  radius,  a,  and  the  stick  zone  radius,  c,  as  derived  by  Hertz 
[17]  and  Mindlin  [19]  respectively  are  given  in  Chapter  2  of  this  work.  The  analytical  values 
for  these  parameters  have  been  calculated  for  each  test  conducted  on  A1  7075-T6  and  Ti-6A1- 
4V  and  will  be  compared  with  the  actual  specimen  values  when  possible.  The  results  of  the 
comparison  are  given  in  the  next  section. 

A  second  objective  of  this  observational  work  was  to  determine  the  location  of  crack  initia¬ 
tion.  To  achieve  this  objective  the  spherical  fretting  pads  and  the  unbroken  fatigue  specimens 
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Figure  3-5:  Plot  of  effective  stress  due  to  contact  versus  number  of  cycles  to  failure  for  Ti-6A1-4V 
(ab  ~  300  MPa). 

have  been  analyzed.  The  fretting  pads  provide  an  opportunity  to  see  the  influence  of  the  con¬ 
tact  stresses  alone  on  crack  initiation,  due  to  the  fact  that  they  are  not  subjected  to  an  applied 
bulk  stress. 

3.4.1  A1  7075-T6 

Not  all  of  the  specimens  from  the  33  A1  7075-T6  tests  shown  in  Table  3.1  were  available  to  be 
analyzed,  only  about  one-half  of  the  tests  had  pads  or  specimens  that  could  be  analyzed.  The 
fretting  pads  provided  nearly  all  of  the  observations  because  when  the  specimens  failed  they 
usually  failed  at  the  fretting  scar.  It  was  very  difficult  to  determine  the  fretting  parameters  a, 
and  c/a  for  many  of  the  pads  and  in  some  cases  it  could  not  be  done.  The  stick  and  slip  zones 
could  not  be  distinguished  for  many  of  the  pads  that  correspond  to  specimen  failure.  Table 
3.9  shows  the  comparison  between  the  predicted  (Pred.  in  the  table)  and  actual  (Act.  in  the 
table)  values  for  the  fretting  geometry,  for  those  specimens  where  the  comparison  was  possible. 
L  and  R  correspond  to  the  left  and  right  fretting  pads  respectively.  Test  33,  which  did  not 
produce  specimen  failure  also  displays  the  values  for  the  specimen  and  are  denoted  by  “(spec)” 
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Figure  3-6:  Picture  of  the  left  fretting  pad  used  for  test  33  with  the  approximate  size  of  the 
stick  zone  shown  (A1  7075-T6). 

in  the  table.  It  can  be  seen  from  this  table  that  the  actual  contact  radius  is  in  general  larger 
than  that  predicted  by  Equation  2.1.  In  some  cases  the  difference  is  quite  substantial.  The 
actual  c/a  ratios  are  fairly  close  to  the  predicted  values,  indicating  that  the  value  used  for  ji  is 
accurate. 

As  discussed  earlier,  a  scanning  laser  microscope  was  used  to  observe  these  specimens  and 
to  produce  images  of  the  scars.  Figure  3-6  shows  the  fretting  scar  on  the  left  pad  used  for  test 
33.  The  stick  and  slip  regions  are  clearly  evident  and  the  eccentricity  caused  by  the  applied 
bulk  can  also  be  seen.  The  corresponding  fretting  fatigue  specimen  is  shown  in  Figure  3-7.  The 
fretting  scar  on  the  specimen  is  not  as  visible  as  the  scar  on  the  pad.  The  fretting  geometry  is 
not  as  easy  to  distinguish  for  specimens  that  were  cycled  to  failure.  It  is  difficult  to  determine 
even  the  contact  radius  of  the  pads  for  the  failed  specimens  and  nearly  impossible  to  determine 
the  stick  zone  radius. 

The  other  objective  of  this  section  was  to  determine  the  location  of  cracking.  No  cracking 
was  evident  in  Figures  3-6  and  3-7.  It  would  have  been  extremely  beneficial  if  cracking  could 
have  been  observed  in  Figure  3-6,  because  the  stick  and  slip  zones  are  readily  identifiable.  As 
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Test 

c/a  (Pred.) 

c/a  (Act.) 

Cracking 

5L> 

207 

200 

0.73  . 

* 

no 

5R 

207 

186 

0.73 

* 

yes 

12R 

159 

247 

0.58 

* 

yes 

15L 

213 

325 

0.55 

* 

no 

15R 

213 

305 

0.55 

* 

no 

17L 

238 

248 

0.74 

* 

no 

17R 

238 

264 

0.74 

* 

no 

18L 

177 

210 

0.72 

* 

no 

19L 

213 

264 

0.72 

* 

no 

19R 

213 

249 

0.72 

* 

no 

21L 

213 

226 

0.72 

* 

no 

21R 

213 

202 

0.72 

* 

no 

22L 

213 

240 

0.72 

0.64 

no 

22R 

213 

241 

0.72 

0.65 

no 

23L 

213 

230 

0.72 

0.51 

no 

23R 

213 

225 

0.72 

0.67 

no 

28L 

213 

362 

0.72 

0.79 

no 

28R 

213 

* 

0.72 

* 

no 

29  L 

213 

* 

0.72 

* 

indefinite 

29  R 

213 

* 

0.72 

* 

indefinite 

30  L 

213 

* 

0.72 

* 

yes 

30  R 

213 

* 

0.72 

* 

yes 

31  L 

213 

241 

0.72 

0.65 

yes 

31  R 

213 

265 

0.72 

* 

yes 

32  L 

213 

239 

0.72 

* 

yes 

32  R 

213 

279 

0.72 

* 

yes 

33  L 

213 

238 

0.72 

0.61 

no 

33  R 

213 

257 

0.72 

0.61 

no 

33  L  (spec) 

213 

253 

0.72 

0.69 

no 

33  R  (spec) 

213 

244 

0.72 

0.78 

no 

Table  3.9:  Comparison  of  the  analytical  versus  actual  values  for  the  fretting  geometry,  A1 
7075-T6  (*  indicates  value  could  not  be  determined). 


Figure  3-7:  Picture  of  the  left  side  of  the  specimen  used  for  test  33  with  the  approximate  size 
of  the  stick  zone  shown  (A1  7075-T6). 

Table  3.9  shows,  cracking  was  found  on  eight  of  the  fretting  pads.  The  cracking  found  on  the 
left  fretting  pad  for  test  32  is  shown  in  Figure  3-8.  The  cracking  is  on  the  right  side  of  the 
picture  and  appears  to  have  started  approximately  at  the  edge  of  contact.  The  stick  and  slip 
zones  cannot  be  distinguished,  so  the  exact  location  of  crack  initiation  cannot  be  determined 
from  this  figure.  It  is  interesting  that  cracking  was  found  on  the  fretting  pads  used  for  test 
32,  but  no  evidence  of  cracking  can  be  found  on  the  pads  or  specimen  used  for  test  33,  The 
two  tests  had  exactly  the  same  contact  loading  and  the  difference  in  the  applied  bulk  stress  was 
only  7  MPa.  It  would  seem  logical  that  fretting  would  have  initiated  cracking  in  both  tests 
and  an  arrested  crack  would  be  produced  in  test  33  due  to  the  lower  bulk  stress,  but  this  was 
not  found. 

3.4.2  Ti-6A1-4V 

The  results  of  the  observations  performed  on  the  Ti-6A1-4V  test  specimens  are  shown  in  Table 
3.10.  Again  it  was  very  difficult  to  observe  the  fretting  geometry  for  the  specimens  which 
were  cycled  to  failure;  however,  test  5  provided  an  interesting  specimen.  Even  though  the 
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Figure  3-8:  Picture  of  cracking  on  the  left  fretting  pad  used  for  test  32  with  the  approximate 
circular  area  of  contact  shown  (A1  7075-T6). 

specimen  was  cycled  to  failure,  the  fracture  surface  went  above  the  scar  on  the  left  side  of  the 
specimen  and  preserved  the  scar.  The  ratio  of  c/a  was  calculated  using  a  coefficient  of  friction, 
^  =  0.5,  which  produced  fairly  good  agreement  between  the  analytical  and  calculated  values 
of  c/a.  The  predicted  values  for  the  radius  of  contact  were  also  in  reasonably  good  agreement 
with  the  actual  values  observed  for  the  tests. 

The  scanning  laser  microscope  was  again  used  to  produce  photos  of  the  pads  and  specimens. 
The  left  fretting  pad  used  in  test  5  is  shown  in  Figure  3-9.  The  slip  and  stick  zones  are  clearly 
visible  in  this  figure,  which  is  unusual  for  a  test  that  went  to  specimen  failure.  The  left  side 
of  the  specimen  used  in  test  5  is  shown  in  Figure  3-10.  This  is  the  specimen  discussed  earlier 
where  the  fracture  plane  was  above  the  fretting  scar  on  the  left  side  of  the  specimen. 

It  is  interesting  that  cracking  cannot  be  observed  in  Figure  3-10.  It  seems  reasonable 
that  cracking  would  initiate  on  both  sides  of  the  specimen  since  both  sides  are  under  the  same 
loading  conditions.  Cracking  w*as  observed  in  the  fretting  pads  of  test  2.  The  cracking  in  the 
right  fretting  pad  used  in  test  2  is  shown  in  Figure  3-11.  Cracking  can  be  seen  on  the  left  edge 
of  the  scar.  Cracking  was  also  found  on  the  opposite  side  of  the  scar,  but  it  cannot  be  seen  in 


Test 

a  (n m>  Pred.) 

a  (fx m,  Act.) 

c/a  (Pred) 

c/a  (Act.) 

Cracking 

1  L 

195 

200 

0 

* 

no 

1R 

195 

241 

0 

0.34 

no 

2  L 

2  R 

195 

195 

216 

247 

0.39 

0.39 

0.47 

* 

3  L 

195 

280 

0.75 

* 

no 

3  R 

195 

* 

0.75 

* 

no 

5  L 

195 

195 

0.78 

0.76 

no 

5  R 

195 

230 

0.78 

* 

no 

5  L  (spec) 

195 

196 

0.78 

0.77 

no 

5  R  (spec) 

195 

175 

0.78 

* 

no 

6  L 

195 

252 

0.75 

* 

no 

6  R 

195 

243 

0.75 

* 

no 

7  L 

195 

252 

0.43 

* 

no 

7  R 

195 

* 

0.43 

* 

no 

8  L 

164 

189 

0 

* 

no 

8  R 

164 

302 

0 

* 

no 

9  L 

164 

233 

0.74 

0.64 

no 

9  R 

164 

197 

0.74 

0.65 

no 

Table  3.10:  Comparison  of  the  analytical  versus  actual  fretting  geometry  (Ti-6A1-4V). 


Figure  3-10:  Picture  of  the  left  side  of  the  specimen  used  for  test  5  (Ti-6A1-4V). 

this  photo.  The  crack  appears  to  have  begun  near  the  edge  of  contact;  however,  the  slip  and 
stick  zones  are  not  distinguishable  and  the  exact  location  of  crack  initiation  is  not  certain. 
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Figure  3-11:  Picture  of  cracking  in  the  right  fretting  pad  used  for  test  2  (Ti-6A1-4V). 


Chapter  4 


Fretting  Fatigue  Life  Prediction 
Analysis 


4.1  Introduction 

There  are  two  basic  approaches  used  to  predict  component  life  for  plain  fatigue:  the  damage 
tolerant  approach  and  the  total  life  approach.  The  damage  tolerant  approach  assumes  an  initial 
flaw  exists  in  the  material  and  the  life  of  the  structure  is  the  number  of  loading  cycles  it  takes 
to  propagate  an  initial  flaw  to  a  critical  length  where  failure  occurs.  In  this  method  only  crack 
propagation  is  considered.  The  total  life  approach  considers  the  whole  failure  crack  formation 
process,  from  initiation  to  component  failure.  This  work  uses  both  approaches  to  analyze  the 
fretting  fatigue  phenomenon  to  provide  as  much  insight  into  the  problem  as  possible. 

Developing  a  life  prediction  method  to  be  used  in  fretting  fatigue  applications  is  difficult 
because  the  effects  of  fretting  are  most  damaging  during  the  stages  of  fatigue  crack  growth 
which  have  the  least  developed  methodologies:  the  crack  initiation  stage  and  the  short  crack 
propagation  stage.  Methods  do  not  even  exist  for  plain  fatigue  to  determine  the  number  of 
cycles  required  to  initiate  a  crack  and  propagate  it  the  distance  required  to  use  long  crack 
propagation  techniques.  The  varying  stress  fields  introduced  by  fretting  make  this  task  even 
more  difficult. 

The  deleterious  effects  of  fretting  on  fatigue  were  shown  in  the  preceding  chapter.  The 
applied  bulk  stress  needed  to  cause  fatigue  failure  in  materials  under  fretting  conditions  can 
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Figure  4-1:  Comparison  of  the  S-N  curves  for  plain  fatigue  and  fretting  fatigue  for  A1  7075-T6 
(<rcc  =  337  MPa). 

be  much  lower  than  the  endurance  strength  of  the  material.  The  endurance  strength  of  A1 
7075-T6  is  176  MPa;  however,  under  fretting  conditions  failures  were  observed  at  bulk  stresses 
as  low  as  74  MPa.  The  damaging  effects  of  fretting  on  fatigue  life  is  shown  graphically  in 
Figure  4-1  for  A1  7075-T6.  The  constant  contact  loading  series  of  tests  ( P  —  20  N,  Q  =15  N) 
has  been  used  for  the  fretting  fatigue  S-N  curve  shown  in  the  figure.  Therefore,  the  fretting 
S-N  curve  represents  only  one  contact  loading  condition.  The  curve  could  change  significantly 
for  a  different  set  of  contact  load  parameters.  This  will  be  discussed  in  more  detail  later  in 
this  chapter. 

The  fretting  S-N  curve  has  been  shifted  down  considerably  relative  to  the  plain  fatigue  S-N 
curve  in  Figure  4-1.  Therefore,  cracks  are  initiating  and  growing  to  failure  in  fretting  conditions 
where  there  was  no  failure  for  the  plain  fatigue  conditions,  even  at  higher  alternating  stresses, 
which  demonstrates  the  critical  role  of  fretting  in  crack  initiation.  The  fretting  fatigue  curve 
also  exhibits  an  endurance  limit.  The  existence  of  an  “endurance”  type  limit  for  fretting  fatigue 
is  discussed  in  more  detail  in  the  following  section. 


4.2  Total  Life  Approach 


The  total  life  approach  is  further  subdivided  into  two  categories;  the  stress-life  method  and 
the  strain-life  method.  The  stress-life  method  is  used  in  high  cycle  fatigue  (HCF)  applications 
where  the  applied  stresses  are  maintained  in  the  elastic  range.  The  strain-life  method  is  used 
for  low  cycle  fatigue  (LCF)  applications  where  plastic  strains  make  up  a  significant  portion  of 
the  total  strains  in  the  material.  The  fretting  fatigue  phenomenon  is  a  HCF  problem  and  the 
stresses  calculated  in  the  previous  chapter  were  found  to  be  within  the  elastic  region,  so  the 
stress-life  method  has  been  chosen  for  further  investigation.  The  stress-life  approach  is  still 
widely  used  in  industry  today. 

The  stress-life  approach  is  based  upon  the  work  of  Wohler  performed  in  the  1850’s  and 
1860’s,  who  developed  the  S  -  N  diagram.  The  S  -  N  diagram  is  a  plot  of  the  alternating 
stress,  S,  versus  the  number  of  cycles  to  failure,  N.  The  work  of  Wohler  clearly  showed 
a  relationship  between  the  applied  stress  amplitude  and  the  number  of  cycles  to  failure  and 
indicated  an  existence  of  an  endurance  limit  for  some  materials.  The  endurance  limit  is  the 
applied  alternating  stress  below  which  fatigue  failure  will  not  occur.  Some  materials  do  not 
exhibit  an  endurance  limit.  For  these  materials  the  endurance  strength  has  been  defined  as 
the  stress  amplitude  which  produces  specimen  failure  at  5  x  106  cycles,  or  more  often  shown  as 
2 N  =  1  x  107  reversals. 

A  mathematical  expression  for  the  stress-life  relationship  found  by  Wohler  was  not  produced 
until  1910  by  Basquin  [26].  Basquin  noted  that  a  linear  relation  existed  between  the  applied 
stress  amplitude  and  the  number  of  reversals  to  failure  when  plotted  on  a  log-log  scale.  The 
relationship  is  given  by: 

^  =  aa  =  a'}  (2N/)b  (4.1) 

where  ua  is  the  applied  stress  amplitude,  a ^  is  the  fatigue  strength  coefficient,  b  is  the  fatigue 
strength  exponent  and  2 Nf  is  the  number  of  reversals  to  failure.  The  fatigue  strength  coeffi¬ 
cient,  crp  is  approximately  equal  to  the  true  fracture  strength  of  the  material,  a  f.  The  fatigue 
strength  exponent,  i>,  is  generally  between  -0.05  to  -0.12.  This  equation  has  been  used  quite 
extensively  for  many  years.  The  values  for  and  b  are  known  for  many  different  materials 
and  can  be  found  in  many  text  books,  including  references  [24]  and  [22], 


The  Basquin  relationship  given  in  Equation  4.1  can  be  rearranged  to  give 


Nf  =  0.5 


(4.2) 


A  modification  to  the  equation  above  is  proposed  for  fretted  components  such  that  the  number 
of  cycles  to  failure  is  now  given  by: 


(4.3) 


where  F  is  the  fretting  correction  factor,  which  is  used  to  account  for  the  premature  initiation  of 
cracking  in  fretted  components.  To  be  consistent  with  earlier  notation,  cr j,  has  been  substituted 
in  for  cra.  Equation  4.3  is  very  simple  and  if  its  validity  can  be  proved,  it  would  be  a  valuable 
design  tool.  The  critical  parameter  in  the  modified  Basquin  expression  is  the  fretting  correction 
factor,  F.  The  values  for  the  fatigue  strength  exponent,  b  and  the  fatigue  strength  coefficient, 
c'f  are  taken  to  be  the  same  as  for  plain  fatigue.  They  have  not  been  modified  to  account  for 
fretting.  Equation  4.3  can  be  rearranged  such  that  F  is  given  by: 


crb 


(4.4) 


The  experimental  results  given  in  Table  3.1  have  been  used  to  determine  F  for  each  test.  The 
analysis  of  F  is  separated  into  two  categories;  the  constant  applied  bulk  stress  results  and  the 
constant  contact  stress  results.  This  was  done  to  determine  the  effect  of  the  contact  and  bulk 
stresses  on  F.  F  is  assumed  to  be  a  function  only  of  the  bulk  and  contact  stresses.  The 
constant  bulk  stress  results  are  given  in  Table  4.1.  The  values  of  F  for  the  unfailed  specimens 
were  determined  using  the  number  of  cycles  where  the  tests  were  stopped;  therefore,  these 
values  of  F  are  maximum  values  for  the  given  contact  conditions. 

Figure  4-2  is  a  graphical  representation  of  the  information  displayed  in  Table  4.1.  Tliis 
figure  shows  the  general  trend  of  the  data,  as  the  effective  stress  due  to  contact  is  increased, 
the  fretting  correction  factor,  F,  also  increases.  The  plot  includes  the  unfailed  specimens  in 
order  to  display  the  behavior  of  F  at  low  contact  stresses.  The  F  values  for  the  unfailed 


Test 

Nf  (Cycles) 

<rb(MPa) 

<7ec(MPa) 

F 

1 

449,500 

83 

335 

2.820 

4 

361,000 

83 

366 

2.896 

6 

530,000 

83 

309 

2.758 

7 

3,300,000  (No  failure) 

83 

264 

2.189 

8 

10,680,000  (No  failure) 

83 

248 

1.860 

9 

803,000 

83 

292 

2.617 

10 

2,940,000 

83 

266 

2.223 

11 

3,450,000  (No  failure) 

83 

264 

2.160 

12 

616,500 

83 

285 

2.706 

18 

2,608,000  (No  failure) 

83 

280 

2.250 

19 

549,000 

83 

337 

2.746 

Table  4.1:  The  necessary  data  to  determine  F  and  its  relationship  with  the  contact  stresses. 


Test 

Nf  (Cycles) 

o-b(MPa) 

F 

19 

549,000 

83 

2.746 

28 

3,050,000 

75 

2.449 

29 

390,000 

100 

2.379 

30 

258,063 

125 

2.001 

31 

181,118 

150 

1.748 

32 

2,497,226 

74 

2.545 

1 

449,500 

83 

2.820 

Table  4.2:  Fretting  correction  factors  for  different  bulk  stresses,  constant  contact  loads. 

specimens  are  only  approximations,  however  they  do  represent  the  maximum  value  that  the 
fretting  correction  factor  could  be  for  these  tests.  As  can  be  seen  in  Figure  4-2,  the  initial 
slope  of  the  F  curve  is  fairly  flat  then  rises  sharply  before  trailing  off  again  at  the  high  contact 
stress  values.  The  initial  slope  may  even  be  less,  because  the  F  values  for  the  unfailed  specimens 
are  their  maximum  values.  These  values  may  actually  be  much  lower.  It  must  be  stressed  that 
these  results  are  for  one  value  of  <r\»  a  different  bulk  loading  will  yield  a  different  curve. 

The  values  of  F  for  the  constant  contact  loading  tests  are  given  in  Table  4-3.  Only  the 
tests  that  were  cycled  to  failure  were  included  in  this  table.  Table  4.2  has  also  included  test 
number  one.  The  contact  loading  for  this  test  is  nearly  identical  to  the  others,  P  =  20.8  N 
versus  20.0  X  used  in  the  other  tests.  The  plot  of  F  versus  cr&  is  displayed  in  Figure  4-3  . 

From  Figure  4-3,  it  can  be  seen  that  for  a  given  contact  stress,  F  decreases  as  the  bulk 
loading  is  increased  except  for  the  two  low  bulk  load  cases  (tests  28  &32).  These  two  cases 
will  be  discussed  again  later. 


o  (MPa) 

•c 

Figure  4-4:  Plot  of  the  combination  of  bulk  and  contact  loading  above  which  fretting  fatigue 
failure  will  occur  for  A1  7075-T6. 

The  constant  contact  stress  analysis  and  the  constant  bulk  stress  analysis  will  now  be 
combined  to  provide  an  estimate  of  the  general  behavior  of  F  for  various  contact  and  bulk 
loadings.  The  knock-down  factor  idea  discussed  earlier  can  aid  in  this  effort.  Two  thresholds 
have  been  found  for  the  two  series  of  tests.  Given  ah  —  83  MPa,  specimen  failure  will  not  occur 
if  (Tec  <  264  MPa,  and  given  aec  =  337  MPa,  specimen  failure  will  not  occur  if  <  70  MPa. 
These  two  thresholds  can  be  combined  with  the  endurance  strength  of  the  material  to  produce 
the  failure  envelope  curve  shown  in  Figure  4-4.  Specimen  failure  is  predicted  for  combinations 
of  contact  and  bulk  loadings  which  fall  above  the  curve. 

Figure  4-4  can  now  be  used  to  analyze  Figure  4-2.  For  crec  <  264  MPa  fretting  will  not 
cause  specimen  failure  if  <rb  <  83  MPa;  therefore  F  =  1.0  for  these  loading  conditions.  Figure 
4-2  can  now  be  redrawn  as  seen  in  Figure  4-5. 

The  behavior  of  the  curve  in  Figure  4-5  is  very  interesting.  Fretting  has  no  effect  on  the 
fatigue  life  of  the  specimen  for  a  large  range  of  aec.  However,  once  the  threshold  value  of  aec 
is  crossed,  the  impact  on  the  fretting  correction  factor,  F,  is  quite  significant.  Below  the 
threshold,  the  contact  stresses  are  not  high  enough  to  initiate  and  grow  the  crack  to  the  length 
required  for  long  crack  propagation  by  the  bulk  stress.  Once  the  threshold  contact  stress  is 
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Figure  4-5:  Revised  plot  of  F  versus  aec  for  A1  7075-T6,  a b  =  83  MPa. 

reached  cracks  will  grow  enough  to  be  controlled  by  the  bulk  stress,  but  the  rate  at  which  this 
occurs  is  not  greatly  effected  by  how  far  above  the  threshold  the  contact  stresses  are.  The 
steep  portion  of  the  F  curve  appears  to  be  a  very  narrow  region  surrounding  the  threshold. 
The  slope  of  the  upper  portion  of  the  curve  is  relatively  flat. 

The  constant  contact  loading  data  in  Table  4.2  and  Figure  4-4  can  now  be  used  to  estimate 
the  behavior  of  F  for  other  bulk  loading  conditions  as  shown  in  Figure  4-6.  The  threshold 
values  for  ab  =  100  MPa  and  ab  =  125  MPa  were  determined  using  the  failure  envelope  curve. 
The  shape  of  the  curves  were  assumed  to  be  similar  to  that  found  for  ab  =  83  MPa.  It  can  be 
seen  that  the  F  curves  for  the  100  and  125  MPa  bulk  loading  cases  pass  through  the  previously 
determined  values  of  2.38  and  2.0  at  crec  =  337  MPa. 

Tests  28  &  32  discussed  earlier  provide  results  that  support  the  validity  of  the  curves  pre¬ 
sented  above.  The  F  value  found  for  these  two  cases  fell  below  the  curve  for  the  results  plotted 
in  Figure  4-3.  The  threshold  crec  determined  from  Figure  4-4  is  311  MPa  for  ab  =  74  MPa. 
Therefore,  the  contact  stress  of  337  MPa  is  approaching  the  threshold  value  and  could  be  on 
the  steep  portion  of  the  F  curve,  which  would  result  in  an  F  value  lower  than  the  expected 
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Figure  4-6:  Plot  of  the  predicted  behavior  of  F  versus  aec  for  various  bulk  loadings  (A1  7075-T6). 
peak  value. 

For  a  given  set  of  loading  conditions  Figure  4-6  can  be  used  to  determine  F,  which  can  then 
be  substituted  back  into  equation  4.3  to  determine  the  component  life.  However,  Figure  4-6  is 
only  a  rough  approximation  of  the  behavior  of  F  for  various  loading  conditions.  Many  more 
tests  need  to  be  conducted  to  determine  a  better  estimate  for  the  behavior  of  F.  The  threshold 
contact  stress  values  are  the  most  critical  values  because  it  is  there  where  the  correction  factor 
is  changing  significantly  and  hence  where  fretting  fatigue  life  is  changing  the  most. 

The  advantages  of  this  method  are  that  it  is  simple  in  form  and  utilizes  the  familiar  stress- 
life  approach  widely  used  in  industry  today.  However,  to  implement  this  method  will  require 
considerable  experimental  work  for  many  different  materials.  The  method  also  does  not  provide 
information  on  the  location  of  crack  initiation.  The  method  uses  the  contact  stresses  at  the 
edge  of  contact,  but  this  is  not  always  the  location  of  crack  initiation. 

The  analysis  performed  in  this  section  has  produced  some  very  interesting  results  which  are: 

•  A  rough  approximation  for  the  component  failure  envelope  as  a  function  of  the  bulk  and 
contact  stresses  has  been  determined  and  is  shown  in  Figure  4-4. 
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Figure  4-7:  Plot  of  F  versus  aec  for  Ti-6A1-4V  (crj>  ^  300  MPa). 

#  The  slope  of  the  top  of  the  F  curve  is  relatively  flat  indicating  that  the  fatigue  life  is  not 
greatly  affected  by  changes  in  the  effective  stress  due  to  contact  in  this  region. 

•  F  acts  almost  as  a  step  function  making  the  determination  of  the  threshold  stress  values 
of  great  importance. 

The  same  analysis  was  also  performed  for  the  Ti-6A1-4V  tests  recorded  in  the  previous 
chapter.  Tests  5-9  from  Table  3.7  were  used  to  determine  F  as  a  function  of  acc  for  crb  ~  300 
MPa.  The  results  are  shown  in  Figure  4-7.  This  curve  is  based  on  limited  experimental  data, 
but  the  same  type  of  behavior  as  seen  with  the  aluminum  alloy  seems  to  be  present  for  the 
titanium  alloy.  For  the  given  bulk  loading,  F  =  1  until  aec  ~  390  MPa  then  it  quickly  jumps 
to  ~  1.7. 

4.3  Damage  Tolerant  Approach 
4.3.1  A1  7075-T6 

The  second  life  prediction  approach  used  to  analyze  the  fretting  fatigue  problem  is  the  damage 
tolerant  approach.  Using  the  experimental  results  from  the  preceding  chapter,  this  work 


attempts  to  determine  the  relative  magnitudes  of  the  crack  nucleation  lifetime  and  the  long  crack 
propagation  lifetime.  The  crack  nucleation  lifetime  used  in  tliis  work  is  the  life  spent  initiating 
and  growing  the  fatigue  crack  to  the  length  where  the  fretting  stresses  are  approximately  zero. 
This  approach  is  used  by  Szolwinski  and  Farris  [30]  in  their  fretting  fatigue  life  prediction 
method.  The  long  crack  propagation  life  then  is  only  influenced  by  the  bulk  stress.  Linear 
elastic  fracture  mechanics  can  be  used  to  determine  this  life.  Tliis  long  crack  propagation 
life  can  then  be  compared  with  the  total  fatigue  life  from  the  experiments  to  determine  the 
relative  magnitudes  of  the  two  stages  of  fatigue  life.  The  approach  is  demonstrated  using  the 
A1  7075-T6  experiments  and  is  applied  in  the  following  section  to  the  Ti-6A1-4V  experiments. 

Paris,  Gomez  and  Anderson  [31]  in  1961  and  Paris  and  Erdogan[32]  in  1963  demonstrated 
that  a  power  law  relationship  exists  between  the  crack  growth  increment  da/d  A7  and  the  stress 
intensity  factor  range  A K.  This  relationship  is  given  by: 

=  C(AK)m  (4.5) 

where  C  and  m  are  scaling  constants,  and 

A  A  —  ATmax  ATmjn,  for  Amjn  >  0 
A K  =  ATmax,  for  ATmill  <  0. 

Observations  were  made  of  the  failed  test  specimens  to  determine  the  shape  of  the  propagating 
crack.  It  was  determined  that  the  cracks  could  be  best  characterized  as  surface  elliptical  cracks 
with  a  ratio  of  a/c  —  0.9.  A  diagram  of  this  geometry  is  shown  in  Figure  4-8,  with  the  thickness 
and  width  of  the  fatigue  specimens  also  given. 

The  stress  intensity  factor  equation  for  the  surface  elliptical  crack  geometry  was  determined 
by  Newman  and  Raju  [33]  in  1981.  For  a  finite  plate  subjected  to  a  remote  uniform  tensile 
stress,  S',  the  stress  intensity  factor  is  given  by: 


Figure  4-8:  Diagram  showing  the  elliptical  crack  geometry. 


where 

/a\i.65  a 

<2=  1  +  1.464  (-J  for  -  <  1. 

The  expression  for  F  was  obtained  using  a  systematic  curve- fitting  procedure  [33].  Note  this 
F  is  not  the  same  F  described  in  the  previous  section  and  a  in  these  expressions  refers  to  the 
crack  geometry  and  not  the  contact  radius.  The  expression  for  F  is  given  by: 


where 


F=  MX  +  M2  (I)'  +  MS  (j)4]  Ugfw 

(4.7) 

Mi  =  1.13-0.09  0 

(4.8) 

0.89 

Mi  =  -0.54  +  — — - — — 

*  0.2  +  ( a/c ) 

(4.9) 

1.0  /  a\ 24 

M3  =  °-5-0.65  +  („/c)+14(1-°-e) 

(4.10) 

g  =  1  +  0.1  +  0.35  fj'j  1(1-  sin0)2 

(4.11) 

Parameter 

Value 

m 

4 

c 

8.6x10“ 12 

ao 

850  /im 

a, 

2.5  mm 

t 

2.5  mm 

Table  4.3:  Values  used  in  the  Paris  life  calculations  for  A1  7075-T6.  Consistent  units  for  C  are 
meters  and  MPa. 


(4.12) 

(4.13) 


The  function,  fw  is  the  finite  width  correction. 

The  above  equations  can  now'  be  used  to  find  K  given  any  a.  The  expression  for  the  stress 
intensity  factor  is  now'  substituted  back  into  the  Paris  form,  Equation  4.5  and  the  equation  is 
numerically  integrated  from  an  initial  crack  length  ao  to  the  critical  crack  length  aj  to  find  the 
crack  propagation  life.  Using  0  =  90° ,  the  propagation  life,  Nfp  is  given  by: 


NfP  — 


jl _  ra* 

icr)A  7T2  Ja o 


da 


C(Acr)  7T2  Jao  (0.67)4a2  [l.05  + 0.27(f)2 -0.145(f)4]  [sec  (1#av/f )]‘ 


(4.14) 


The  values  used  in  the  propagation  life  calculations  are  displayed  in  Table  4.3  [35].  The 
value  used  for  ao  is  approximately  four  times  the  contact  radius  for  the  tests  conducted  on  the 
aluminum  alloy.  The  stresses  due  to  contact  decay  away  very  quickly  and  are  approximately 
zero  at  a  depth  of  1.5a.  However,  if  1.5a  is  used  for  ao,  the  calculated  propagation  life  is  greater 
than  the  experimental  life,  W'hich  cannot  be  true.  Therefore,  ao  was  set  equal  to  4a.  The 
value  used  for  ao  must  also  be  long  enough,  such  that  it  includes  at  least  a  few  grain  diameters. 
The  average  grain  size  for  the  A1  7075-T6  specimens  w'as  60  /zm  [23],  so  the  ao  used  in  these 
calculations  encompasses  several  grain  diameters.  As  can  be  seen  from  Table  4.3,  only  one-half 
the  thickness  of  the  specimen  was  used  in  the  calculations.  This  is  due  to  the  presence  of  two 
cracks  in  the  specimens.  When  the  fracture  surfaces  of  the  failed  specimens  were  analyzed, 
two  cracks  were  observed.  One  on  each  side  of  the  fatigue  specimen. 

Using  Equation  4.14,  the  long  crack  propagation  life  wras  calculated  for  each  test,  using  ab 


for  A  a.  The  stress  intensity  factor  range  was  also  calculated  for  each  test  using  Equation  4.6, 
at  a  =  o0,  and  S  =  crb  to  examine  if  AK  >  AK0 ,  the  threshold  stress  intensity  factor  range 
at  this  crack  length.  The  A K0  for  A1  7075-T6  at  R  =  -1  appears  to  be  approximately  2-3 
MPa^/ni-  This  value  was  determined  from  a  plot  in  [34].  The  plot  did  not  contain  R  =  —  1 
crack  growth  data;  however,  using  the  displayed  trends  for  the  other  R  ratios  an  estimate  for 
the  AAo  for  R  =  -1  was  made.  The  calculated  values  of  the  propagation  life,  NfP  and  the 
stress  intensity  factor  range  at  a  =  ao  for  each  test  are  given  in  Table  4.4. 

As  can  be  seen  from  Table  4.4,  the  long  crack  propagation  life  is  a  very  significant  portion 
of  the  fretting  fatigue  life  for  many  of  the  failed  specimens.  For  these  specimens,  a  very  small 
portion  of  the  life  was  used  to  nucleate  cracks.  The  calculated  A K  values  in  Table  4.4  provide 
insight  into  some  of  the  tests  where  component  failure  did  not  occur.  It  can  be  seen  that  for 
Tests  14-17  and  Tests  21-24,  A K  <  2.5  MPa^/m  for  each  of  these  tests,  which  is  very  close  to 
the  threshold  value  for  this  material. 

4.3.2  Ti-6A1-4V 

The  method  used  for  A1  7075-T6  above  is  applied  to  Ti-6A1-4V  in  this  section.  The  crack 
geometry  for  the  titanium  alloy  specimens  was  very  similar  to  that  found  above.  Again  elliptical 
surface  cracks  were  found;  however,  cracking  was  found  on  only  one  side  of  the  specimen  for  this 
material.  The  values  used  to  calculate  the  propagation  life  for  this  material  are  shown  in  Table 
4,5  (a/c  corresponds  to  the  elliptic  crack  geometry).  The  values  used  for  ao  are  approximately 
1.5a  and  4a.  The  first  value  of  ao  ^  1.5a,  corresponds  to  the  depth  at  which  the  tensile  stresses 
due  to  contact  are  approximately  zero.  The  second  value  of  ao  ~  4a  is  used  to  be  consistent 
with  the  calculations  performed  earlier  for  A1  7075-T6. 

The  values  for  m  and  C  were  determined  using  curves  found  in  [37]  and  [36] .  The  fatigue 
crack  growth  curve  in  [36]  was  for  solution  treated  and  over-aged  (STOA),  Bi- modal  Ti-6A1-4V, 
while  that  found  in  [37]  was  for  Bi-modal  Ti-6A1-4V  with  no  heat  treatment.  Both  curves  were 
made  using  R  =  0.1.  The  values  found  for  C  and  m  for  the  two  cases  are  shown  in  Table  4.6. 
The  experimental  results  displayed  in  Chapter  3  are  for  Bi-modal  Ti-6A1-4V  /with  no  special 
heat  treatment,  at  R  =  —  1.  The  value  of  m  for  this  section  was  set  to  5.0,  since  the  applied 
stress  ratio  has  a  small  effect  on  m.  The  value  of  C  was  determined  by  approximating  the 


Table  4.4:  Long  crack  propagation  life  for  A1  7075-T6. 


Parameter 

Value 

a/c 

1.0 

t 

5.0  mm 

2  W 

5.0  mm 

m 

5.0 

C 

7.0  xl0~14 

ao 

290  jam,  780  /xm 

°/ 

5  mm 

Table  4.5:  Values  used  for  the  propagation  life  calculations  for  Ti-6A1-4V. 


STOA,  Bi-modal  Ti-6A1-4V 

Bi-modal  Ti-6A1-4V 

C 

8.13x10“ 14 

m 

5.12 

Table  4.6:  Values  of  C  and  m  for  STOA,  Bi-modal  Ti-6A1-4V  and  Bi-modal  T1-6A1-4V  (units 
are  meters  and  MPa). 

effect  of  R  on  C. 

Using  the  values  displayed  in  Table  4.5,  the  propagation  life  was  calculated  for  the  nine 
tests  shown  in  Table  3.7.  The  propagation  lives  corresponding  to  the  two  different  values  of 
a0  are  shown  in  Table  4.7.  As  can  be  seen  from  the  table,  the  propagation  life  correlates  very 
well  with  the  actual  experimental  life  of  the  failed  specimens  when  ao  =  1.5a  is  used.  Using 
ao  =  1.5a,  the  difference  in  fatigue  life  for  the  various  tests  can  be  attributed  almost  entirely  to 
the  difference  in  propagation  life.  Using  ao  =  4a  results  in  much  lower  values  of  Nfp.  Table 
3.7  also  displays  the  calculated  AK  values  at  ao  =  1.5a.  The  threshold  stress  intensity  factor 
range,  AKq,  is  approximately  3-7  MPa^/m  [37]. 


Test 

N /  (Experiment) 

N/p(ao  =  1.5a) 

N /p(ao  =  4.0a) 

AK  (MPay/m,ao  =  l-5a) 

1 

691,612 

1,010,000 

4.79 

2 

109,574 

■ 

| 

8.51 

3 

107,879 

8,130 

8.38 

4 

31,540  (interrupted) 

100,000 

13,300 

7.60 

5 

304,456 

247,000 

32,900 

6.34 

6 

474,656 

239,000 

31,800 

6.38 

7 

143,683 

156,000 

20,700 

6.96 

8 

258,576 

239,000 

31,800 

6.38 

9 

3,257,000(no  failure) 

239,000 

31,800 

6.38 

Table  4.7:  Ti-6A1-4V  long  crack  propagation  life. 


Chapter  5 


Palliatives 


5.1  Introduction 

The  work  performed  in  the  previous  chapters  provide  the  necessary  background  for  a  discussion 
of  fretting  fatigue  palliatives.  The  methods  developed  earlier  provide  the  criteria  that  can  be 
used  to  determine  how  and  why  certain  palliatives  work  while  others  do  not. 

This  work  will  focus  on  two  fretting  fatigue  palliatives:  shot  peening  and  coatings.  Both 
have  been  shown  to  be  beneficial  in  increasing  fretting  fatigue  life.  Chivers  and  Gordelier  [41] 
provide  a  detailed  study  of  several  possible  palliatives  including  shot  peening  and  coatings. 

5.2  Shot  Peening 

Shot  peening  has  been  used  for  many  years  as  a  way  to  improve  the  fatigue  properties  of 
components  in  plain  fatigue.  The  benefits  of  shot  peening  were  found  quite  by  accident  by 
the  Buick  Motor  Division  of  General  Motors  Corporation  in  the  late  1920’s  [42].  It  was  found 
that  valve  springs  that  had  been  blasted  with  grit  to  remove  scale  had  greatly  improved  fatigue 
properties.  It  was  first  thought  that  the  benefits  derived  from  shot  peening  were  due  to  the 
work  hardening  of  the  surface.  However,  it  was  soon  discovered  that  the  benefits  were  caused 
by  the  residual  stresses  introduced  at  the  surface  due  to  shot  peening. 

In  the  shot  peening  process,  relatively  hard  particles  are  shot  against  the  metal  being  treated. 
The  particles  are  generally  spherical  in  shape  and  are  given  enough  velocity  such  that  the  surface 


Figure  5-1:  Plot  showing  approximate  relationship  between  shot  peened  and  unpeened  speci¬ 
mens  for  plain  and  fretting  fatigue. 

of  the  metal  is  indented  during  impact.  The  indentations  occur  due  to  local  plastic  yielding  and 
as  the  surface  material  attempts  to  expand,  it  is  restrained  by  surrounding  material  below  the 
surface  which  was  not  deformed  plastically.  This  restraint  of  the  plastically  deformed  surface 
layer  causes  it  to  be  compressively  strained  and  is  residually  stressed  in  compression. 

There  have  been  numerous  studies  conducted  showing  the  benefits  of  shot  peening  for 
plain  fatigue  conditions.  Shot  peening  has  also  been  demonstrated  to  be  effective  in  fretting 
fatigue  conditions  as  well,  which  is  discussed  in  [38],  [39]  and  [41].  In  many  materials,  shot 
peening  has  been  shown  to  increase  the  fatigue  strength  such  that  the  fretting  fatigue  strength 
is  approximately  the  same  as  the  plain  fatigue  strength.  A  simple  plot  displaying  the  typical 
relationship  between  shot  peened  and  unpeened  specimens  for  fretting  and  plain  fatigue  is 
shown  in  Figure  5-1.  This  figure  is  based  upon  the  work  of  Bignonnet  [39],  which  studied  the 
effects  of  shot  peening  on  fretting  fatigue.  However,  this  figure  is  not  true  for  all  materials.  It 
is  only  meant  to  display  the  benefits  typically  observed  in  shot  peened  materials. 

As  discussed  earlier,  the  benefits  of  shot  peening  are  derived  from  the  residual  compressive 
stresses  introduced  into  the  material  during  the  process.  The  maximum  value  of  the  residual 
compressive  stress,  often  called  the  magnitude,  is  not  affected  by  the  shot  peening  process  as 
long  as  the  shot  used  is  at  least  as  hard  or  harder  than  the  material  being  treated  [43].  The 
magnitude  of  the  residual  compressive  stress  is  material  dependent,  but  has  a  value  of  at  least 
one-half  the  yield  strength  for  that  material  [43].  The  depth  of  the  compressive  layer  is  the 


other  important  parameter  in  shot  peening.  The  depth  of  the  compressive  layer  is  a  function 
of  the  hardness  and  ductility  of  the  peened  material,  the  state  of  strain  in  the  material  when 
the  peening  was  performed  and  the  characteristics  of  the  shot  and  shot  stream  [42].  A  typical 
residual  stress  profile  due  to  shot  peening  and  the  residual  stress  profile  used  in  this  analysis  are 
shown  in  Figure  5-2.  The  residual  stress  profile  used  in  this  analysis  assumes  a  flat  maximum 
compressive  stress,  which  goes  directly  to  zero  at  a  given  depth.  The  assumed  residual  stress 
profile  should  actually  have  a  tensile  region  to  balance  the  residual  compressive  stresses  to 
satisfy  equilibrium.  This  is  shown  in  Figure  5-3,  where  arc  is  the  residual  compressive  stress 
and  art  is  the  residual  tensile  stress  required  for  equilibrium.  For  equilibrium  to  be  satisfied 
area  A\  must  equal  area  A2.  Therefore  the  area  balance  is  given  by: 

(pd)(arc)  =  ( t -pd)(art),  (5.1) 

and  the  ratio  of  the  residual  stresses  is  given  by: 

—  =  77~?T  •  <5'2> 

arc  {t  -  pd ) 

If  (t  —  pd)  (pd)  then  art  — *  0.  It  is  assumed  in  this  work  that  (t  —  pd)  (pd),  therefore  the 
residual  tensile  stress  has  been  neglected.  Note  that  the  actual  maximum  compressive  stress 
occurs  at  some  depth  below  the  surface.  The  actual  depth  of  residual  compressive  stresses  for 
various  materials  is  given  in  [43]  as  w'ell  as  some  residual  stress  profiles. 

The  benefits  derived  from  shot  peening  can  also  become  detrimental  to  the  material  if  the 
residual  compressive  stresses  are  somehow'  relieved.  The  process  which  caused  the  residual 
stresses  to  be  induced  into  the  material  also  damages  the  surface  of  the  material  and  may  even 
introduce  microcracks.  These  microcracks  could  drastically  reduce  fatigue  life  if  the  residual 
compressive  stresses  are  relieved.  This  has  been  showm  to  be  the  case  in  both  plain  fatigue 
and  fretting  fatigue.  Gray,  Wagner  and  Ltitjering  [44]  show  the  detrimental  effects  on  plain 
fatigue  life  of  shot  peening  wdien  the  residual  compressive  stresses  are  thermally  relieved.  The 
fatigue  strength  of  the  material  is  significantly  reduced  for  the  shot  peened,  stress  relieved 
specimens  as  compared  to  the  original  electrolytically  polished  specimens.  Figure  5-4,  from 
[44],  demonstrates  this  effect.  In  Figure  5-4,  SP  is  shot  peened,  EP  is  electrolytically  polished 
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Figure  5-2:  Typical  stress  profile  after  shot  peening  and  the  assumed  stress  profile  used  in  this 
work. 


1  _ AKs-i_ _ , _ 4 

pd  -  peening  depth  _ 


T 


A2 


t  -  material  thickness 


arc  Compressive  Stress  0  °rt  Tensile  Stress 


Figure  5-3:  Residual  stress  equilibrium  diagram. 


Figure  5-4:  Diagram  showing  the  detrimental  effects  of  shot  peening  if  the  residual  stresses  are 
relieved  [44]. 

and  SR  is  stress  relieved.  From  this  figure  it  can  be  seen  that  the  fatigue  strength  for  the 
shot  peened  and  stress  relieved  specimens  is  about  one-half  that  of  the  electrolytically  polished 
specimens  at  lxlO7  cycles,  the  reduction  in  fatigue  strength  is  quite  significant  and  certainly 
cannot  be  ignored.  The  shot  peened  and  stress  relieved  curve  shifts  down  relative  to  the 
electrolytically  polished  specimen  curve  indicating  that  the  difference  between  the  two  curves 
is  in  the  crack  initiation  stage,  just  as  was  found  for  fretting.  The  material  used  in  [44]  was 
Ti-6A1-4V  with  a  fine  lamellar  microstructure  and  the  temperature  was  20  °C  for  this  set  of 
test  results.  The  residual  stresses  were  relieved  thermally  using  a  heat  treatment  of  600  °C  for 
one  hour.  These  results  were  for  plain  fatigue  experiments.  Bignonnet  [39]  displays  similar 
results  for  fretted  specimens.  The  fretting  fatigue  strength  of  the  shot  peened,  stress  relieved 
specimens  was  similar  to  that  of  the  unpeened  specimens.  Waterhouse  [38]  found  similar  results 
in  his  experimental  work.  These  works  demonstrate  the  importance  of  the  compressive  residual 
stresses  and  the  importance  of  not  allowing  them  to  be  relieved. 

Laser  shot  peening  does  not  have  the  shortcomings  discussed  above  for  the  typical  shot 
peening  process.  In  the  laser  shot  peening  process  the  surface  is  unchanged  by  the  peening 
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process.  Shock  waves  are  introduced  into  the  material  that  are  capable  of  plastically  yielding 
the  material  and  producing  compressive  residual  stresses. 

The  residual  stresses  can  be  relieved  either  mechanically  or  thermally.  The  approximate 
temperature  at  which  the  residual  compressive  stresses  begin  to  dissipate  for  various  materials 
can  be  found  in  [43].  The  residual  stresses  are  relieved  mechanically  if  the  total  stresses  in  the 
material  exceed  the  dynamic  elastic  limit  of  the  material.  The  total  compressive  stress  is  made 
up  of  the  shot  peening  residual  stress  and  the  stress  due  to  external  loading.  The  residual 
stress  is  relieved  by  the  amount  that  the  total  compressive  stress  exceeds  the  dynamic  elastic 
limit  of  the  material. 

The  objective  of  this  section  is  find  the  optimum  peening  depth  for  various  fretting  fa¬ 
tigue  loading  conditions.  The  magnitude  of  the  residual  compressive  stresses  are  fixed  by  the 
material,  therefore  the  designer  has  no  control  over  this  aspect  of  the  shot  peening  process. 
However,  the  depth  of  residual  stress  can  be  influenced  by  variations  in  peening  parameters  as 
shown  in  [43].  Therefore,  the  depth  of  peening  becomes  the  variable  design  parameter  that 
can  be  optimized  by  the  designers. 

5.2.1  Material  Yielding  Analysis 

A  systematic  analysis  of  the  shot  peening  process  applied  to  the  sphere  on  flat  surface  contact 
geometry  has  been  conducted  in  this  work.  The  material  used  for  this  analysis  is  Ti-6A1- 
4V.  The  pertinent  material  properties  for  Ti-6A1-4V  are  listed  in  Table  3.6.  The  first  task 
was  to  find  the  conditions  W'hen  the  residual  compressive  stresses  would  begin  to  be  relieved. 
According  to  [43],  the  compressive  stresses  begin  to  dissipate  thermally  at  480  C  for  Ti-6A1-4V. 
Mechanical  stress  relieving  begins  to  occur  when  the  material  stresses  exceeds  the  elastic  limit 
of  the  material.  The  stress  fields  for  a  sphere  fretting  against  a  flat  surface  have  already  been 
discussed  in  Chapter  2  of  this  work.  The  surface  stresses  for  this  geometry  for  A17075-T6  are 
displayed  in  Figure  2-2.  The  stress  profile  for  Ti-6A1-4V  is  almost  identical,  the  magnitudes 
are  different,  but  the  shape  of  the  curves  are  the  same. 

The  analysis  performed  in  this  chapter  assumes  the  onset  of  global  sliding.  This  condition 
provides  the  highest  stresses  possible  for  a  given  contact  loading.  If  global  sliding  is  assumed, 
the  normal  load,  P,  becomes  the  critical  contact  parameter  since  the  tangential  load,  Q ,  will 


always  be  equal  to  fiP.  The  loading  which  causes  the  total  effective  stress,  (Tetoti  to  become 
equal  to  the  yield  strength  is  when  material  yielding  is  assumed  to  occur  and  is  shown  in  Figure 
5-5.  The  data  used  to  make  this  figure  were  determined  using  a  MATLAB  script  file.  A 
MATLAB  routine  was  written  which  first  scans  the  surface  of  the  material,  then  incrementally 
scans  below  the  contact  surface  calculating  the  stresses  at  each  increment,  producing  a  matrix 
of  stresses.  The  surface  is  scanned  along  the  y  =  0  line,  where  maximum  material  stresses  are. 
The  increment  used  in  the  x-direction  is  0.02a.  The  increment  in  the  z-direction  is  0.025a. 
The  program  then  finds  the  maximum  effective  stress  from  the  matrix  of  stresses  determined 
during  the  scan.  The  stresses  at  the  contact  edges  are  calculated  separately.  The  maximum 
stress  is  then  compared  to  the  yield  stress.  If  the  maximum  effective  stress  is  less  than  the 
yield  stress,  the  contact  loading  is  increased  and  the  stress  calculations  are  performed  again. 
The  process  is  continued  until  the  contact  loading  which  produces  yielding  is  determined.  This 
was  done  for  five  different  coefficients  of  friction  and  eight  different  bulk  stress  conditions.  The 
stresses  at  the  edge  of  contact  are  calculated  separately  because  the  stresses  are  indeterminate 
at  the  contact  origin  (x  =  y  =  z  =  0)  and  the  origin  would  be  included  in  the  surface  scan  if 
both  contact  edges  are  also  included.  The  program  used  for  this  analysis  is  shown  in  Appendix 
A. 

In  Figure  5-5  the  yielding  curves  have  been  plotted,  which  show  for  a  given  bulk  stress, 
<jb  and  a  given  coefficient  of  friction,  /z,  the  minimum  contact  pressure,  Po,  that  is  needed  to 
cause  the  total  effective  stress  to  exceed  the  yield  stress  of  the  material.  For  example  if  ji  —  0.7 
and  crb  —  0.1  cry  then  the  maximum  contact  pressure  which  produces  yielding  is  ~0.28 cry.  This 
example  is  shown  in  Figure  5-5.  The  maximum  contact  pressure  is  normalized  by  the  yield 
strength  of  the  material  in  this  plot.  A  residual  stress  of  -0.85cry  is  used  for  this  study.  This 
value  for  the  residual  compressive  stress  was  determined  using  an  actual  shot  peening  stress 
profile  from  [44]  and  a  graph  of  maximum  compressive  residual  stress  versus  tensile  strength  in 

[43]. 

As  can  be  seen  in  Figure  5-5,  the  contact  loading  which  causes  material  yielding  increases 
as  the  applied  bulk  stress  increases  until  ab  —  0.4<ry.  The  leading  edge  of  contact,  where  the 
compressive  stresses  due  to  contact  are  very  high,  is  the  location  of  yielding  for  applied  bulk 
stresses  less  than  one-half  the  yield  strength.  The  location  of  yielding  shifts  to  the  trailing 


Figure  5-5:  The  combination  of  contact  and  bulk  loading  which  causes  the  total  effective  stress 
to  equal  the  yield  stress  plotted  versus  the  coefficient  of  friction. 

edge,  where  the  tensile  stresses  due  to  contact  are  very  high,  when  crb  >  O.oay.  Therefore, 
yielding  is  controlled  by  the  compressive  stresses  until  the  applied  bulk  stress  reaches  a  value 
of  one-half  the  yield  strength.  This  explains  why  the  maximum  contact  pressure  which  causes 
yielding  actually  increases  with  increasing  bulk  stress  when  ab  <  0.4<ry .  Figure  5-5  also  shows 
the  effect  of  friction.  As  the  coefficient  of  friction  increases,  the  contact  pressure  needed  to 
cause  yielding  decreases.  It  can  also  be  seen  that  the  contact  loads  needed  to  cause  yielding  are 
relatively  low  for  the  coefficient  of  friction  typical  for  Ti-6A1-4V  (0.4-0. 5);  therefore,  residual 
stress  relieving  due  to  yielding  cannot  be  ignored.  Yielding  first  occurs  at  the  surface  for  /i  >  0.3. 
Since  yielding  usually  occurs  first  at  the  surface,  the  depth  of  yielding  must  also  be  found. 

5.2.2  Depth  of  Yielding 

The  depth  of  yielding  must  be  found  to  determine  to  what  depth  the  residual  stresses  are  being 
relieved.  The  von  Mises  effective  stress  is  again  used  to  determine  the  yield  conditions,  i.e. 
when  (Tet.ot  >  cry  the  material  is  considered  to  have  yielded.  A  MATLAB  routine  was  also 
written  to  perform  this  task.  The  routine  used  to  perform  this  task  is  similar  to  the  previous 


Figure  5-6:  Plot  of  yielding  depth,  z  versus  maximum  contact  pressure  for  various  bulk  loadings 

0*  =  1.0). 

one.  Again  the  surface  is  scanned  using  an  increment  of  0.02a  in  the  x-direction  with  y  =  0 
for  all  points.  The  stresses  are  then  determined  below  the  surface  using  an  increment  of  0.02a 
in  the  z-direction,  which  results  in  a  matrix  of  stresses.  The  stresses  at  the  edge  of  the  contact 
are  not  calculated  separately  in  this  routine  because  the  stress  profile  is  much  smoother  at  the 
edges  below  the  contact  surface  as  shown  in  Figure  2-3.  The  depth  of  maximum  stress  is  found 
and  this  depth  became  the  starting  point  for  the  iteration  to  find  the  depth  of  yielding.  The 
depth  of  maximum  stress  is  always  the  surface  for  (i  >  0.3.  The  scanning  process  begins  again 
at  the  depth  just  determined  using  a  reduced  increment  in  the  z-direction  of  0.01a.  The  depth 
is  increased  until  the  effective  stress  is  no  longer  greater  than  the  yield  strength  of  the  material. 
Figure  5-6  displays  the  results  of  this  analysis  for  ^  =  1.0.  The  depth,  z  is  normalized  by  the 
radius  of  contact,  a,  the  maximum  contact  pressure,  Po>  is  normalized  by  the  yield  strength  of 
the  material,  ay.  Again  the  residual  stress  used  in  this  analysis  w*as  <7residuai  =  —0.85 cry.  The 
MATLAB  script  file  used  in  this  section  is  shown  in  Appendix  B. 

The  results  displayed  in  Figure  5-6  follow  the  expected  trends,  the  depth  of  yielding  increases 
as  the  contact  load  is  increased  and  decreases  as  the  bulk  loading  is  increased  until  cr&  ^  0.5ay . 
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Figure  5-7:  Depth  of  yielding  versus  maximum  contact  pressure  for  various  bulk  stresses  (jj,  = 
0.8). 

The  depth  of  yielding  decreases  with  increasing  bulk  stress  due  to  the  large  residual  compressive 
stress  assumed  in  the  analysis.  This  same  analysis  was  repeated  for  three  more  coefficients  of 
friction,  fx  =  0.8,  fi  —  0.5  and  fi  =  0.3.  The  results  for  the  three  other  coefficients  of  friction 
are  shown  in  Figures  5-7,  5-8  and  5-9  respectively.  As  can  be  seen  in  these  figures,  the  depth 
of  yielding  is  reduced  significantly  as  the  coefficient  of  friction  is  lowered,  which  again  points 
to  the  possible  benefits  of  coatings  in  fretting  situations.  When  the  coefficient  of  friction  is 
reduced  to  0.3,  there  is  very  little  yielding  even  at  very  high  contact  loadings  for  crb  >  0.4. 

The  depth  of  yielding  curves  developed  in  this  section  can  be  very  useful  for  design,  specif¬ 
ically  in  the  material  treatment  stage  of  design.  Designers  could  use  these  plots  with  the 
known  contact  and  bulk  loadings  and  the  coefficient  of  friction  and  very  quickly  determine 
the  minimum  depth  of  shot  peening  to  ensure  that  the  compressive  residual  stresses  are  not 
mechanically  relieved.  For  example,  if  fi  =  0.8,  Po/&y  =  0.8,  and  ab  =  O.Ioy,  the  depth 
of  yielding,  z/a  ~  0.35.  Therefore,  the  minimum  depth  of  shot  peening  to  avoid  mechanical 
relieving  of  the  residual  stresses  is  z/a  >  0.35.  This  example  is  shown  in  Figure  5-7.  The 
depth  of  yielding  curves  provide  important  design  information,  but  they  only  provide  part  of 


the  information  needed  to  determine  the  optimum  shot  peening  depth.  These  curves  provide 
the  minimum  shot  peening  depth  to  avoid  residual  stress  relieving  but  they  do  not  provide 
information  pertaining  to  the  tensile  stresses  at  the  crack  tip.  An  additional  analysis  is  needed 
to  provide  this  information,  which  is  discussed  in  the  next  section. 

5.2.3  Maximum  Principal  Stress  Analysis 

The  predicted  success  of  shot  peening  is  based  upon  two  basic  assumptions.  First,  a  crack  will 
not  propagate  unless  a  tensile  stress  forces  it  open  near  the  tip  [43].  The  second  assumption 
is  that  a  crack  tip  will  not  open  while  a  compressive  force  is  acting  upon  it  [43].  These 
assumptions  provide  the  criteria  that  is  used  to  determine  the  optimum  shot  peening  depth. 
Using  the  two  assumptions  stated  above,  the  depth  of  peening  to  prevent  crack  propagation 
should  be  the  depth  at  which  the  maximum  principal  stress,  cr i  <  0.  Therefore,  inside  this 
layer  of  material  affected  by  shot  peening,  no  region  of  tensile  stresses  will  exist.  This  does 
not  mean  that  the  entire  depth  of  the  material  will  have  no  tensile  stresses.  This  only  applies 
in  the  shot  peened  region,  but  that  is  enough  to  stop  the  crack  from  propagating. 

A  MATLAB  program  similar  to  the  two  discussed  above  was  written  and  is  shown  in 
Appendix  C.  The  surface  is  scanned  along  the  x-direction  determining  the  stresses  at  each 
point  using  an  increment  in  the  x-direction  of  0.02a,  again  y  —  0  for  all  points.  At  each  point 
the  maximum  principal  stress,  <7i,  is  determined  creating  a  vector  of  maximum  principal  stresses. 
The  maximum  of  this  vector  is  then  found  and  compared  to  zero.  If  the  maximum  of  the  vector 
is  greater  than  zero,  then  the  depth  is  increased  by  0.01a  and  the  process  continued  until  the 
maximum  principal  is  less  than  or  equal  to  zero.  This  depth  then  became  the  required  shot 
peening  depth.  The  shot  peening  depth,  z,  which  produced  a  maximum  principal  stress  less 
than  or  equal  to  zero  was  found  for  many  different  loading  conditions.  Just  as  in  the  yielding 
depth  curves,  the  depth  is  plotted  versus  maximum  contact  pressure  for  several  different  values 
of  applied  bulk  stress.  Figure  5-10  displays  the  depth  z,  at  which  cr\  <  0,  for  fx  =  1.0.  Again 
the  depth  has  been  normalized  by  the  contact  radius,  a  and  the  maximum  contact  pressure, 
Po,  is  normalized  by  the  yield  strength  of  the  material.  This  analysis  was  repeated  for  three 
other  coefficients  of  friction,  (x  =  0.8,  fx  =  0.5  and  (x  —  0.3.  The  curves  for  /i  =  0.8  and  jx  =  0.5 
are  shown  in  Figures  5-11  and  5-12  respectively.  The  curves  for  fx  —  0.3  are  not  displayed 


Figure  5-10:  The  depth  at  which  the  maximum  principle  stress  is  less  than  or  equal  to  zero 
plotted  versus  maximum  contact  pressure,  for  magnitudes  of  bulk  stress  (fi  =  1.0). 

because  the  required  peening  depth  was  less  than  for  0.10a  for  even  the  highest  load  conditions 
analyzed.  The  plots  again  show  the  advantages  of  lowering  the  coefficient  of  friction.  There  is 
a  significant  reduction  in  z/a  as  the  coefficient  of  friction  is  reduced.  As  expected  the  required 
shot  peening  depth  increases  as  the  applied  bulk  stress,  a t  is  increased.  As  seen  in  the  figures, 
the  increase  in  peening  depth  with  increasing  at  becomes  quite  significant  when  at  >  OAay- 
The  plots  generated  in  this  section  cannot  be  used  in  isolation,  they  must  be  integrated  with 
the  plots  from  the  previous  two  sections  of  this  analysis.  Together,  these  sections  provide  the 
designer  the  information  necessary  to  determine  the  optimal  shot  peening  depth  for  the  given 
loading  and  material  conditions. 

5.2.4  Integrated  Yielding  and  Maximum  Principal  Stress  Analysis  (Opti¬ 
mum  Shot  Peening  Depth) 

The  depths  determined  from  the  two  preceding  sections  provide  two  different  threshold  val¬ 
ues.  The  depths  output  from  the  material  yielding  analysis  provide  the  minimum  depth  to 
which  the  specimen  must  be  peened  to  prevent  mechanical  relieving  of  the  residual  compres- 
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Figure  5-11:  The  depth  at  which  the  maximum  principal  stress  is  less  than  or  equal  to  zero 
plotted  versus  maximum  contact  pressure  for  various  magnitudes  of  bulk  stress  (fj,  =  0.8). 


Figure  5-12:  The  depth  at  which  the  maximum  principal  stress  is  less  than  or  equal  to  zero 
plotted  versus  the  maximum  contact  pressure  for  various  magnitudes  of  bulk  stress  (jjl  =  0.5). 


sive  stresses.  The  depth  of  yielding  analysis  determines  the  minimum  peening  depth  to  ensure 
that  the  residual  compressive  stresses  remain  in  the  material.  The  maximum  principal  stress 
analysis  assumes  the  residual  compressive  stresses  are  still  present  and  determines  the  peening 
depth  required  to  prevent  crack  growth.  Therefore,  the  depth  of  yielding  analysis  provides  a 
minimum  shot  peening  depth,  but  does  not  provide  all  information  necessary  to  determine  the 
optimal  peening  depth.  The  remaining  information  is  provided  by  the  maximum  principal 
stress  analysis,  which  determines  the  peening  depths  at  which  o\  <  0.  This  analysis  assumes 
the  crack  will  not  grow  in  a  region  of  pure  compressive  stresses.  The  two  analysis  must  be 
combined  to  form  the  complete  picture,  which  has  been  done  in  Figure  5-13  for  ji  =  1.0  and 
ab  =  0.4cry.  As  can  be  seen  from  Figure  5-13  the  limiting  curve  depends  upon  the  contact 
loading,  a  cross-over  point  occurs  at  (Po/cy)  c*  1.15.  If  (Po/^y)  <  1.15  the  maximum  prin¬ 
cipal  stress  analysis  is  the  limiting  case  and  controls  the  optimum  peening  depth,  however  if 
(Po/cry)  >  1.15  the  depth  of  yielding  analysis  becomes  limiting  and  determines  the  required 
peening  depth.  This  example  shows  the  importance  of  examining  both  criteria.  An  example 
on  how  to  use  the  figure  is  also  displayed  in  Figure  5-13.  If  (Po/cry)  =  0.8  then  the  optimum 
peening  depth,  z/a  ~  0.4.  The  designer  may  choose  to  shot  peen  to  a  slightly  greater  depth 
to  ensure  no  tensile  stresses  exist  in  the  material  at  that  depth. 

A  second  example  is  given  in  Figure  5-14  for  /z  =  0.8  and  crb  —  O.lay.  As  can  be  seen  in 
tliis  figure,  the  intersection  between  the  two  curves  occurs  at  a  much  lower  contact  pressure 
than  in  the  previous  example.  Contrary  to  the  previous  example,  the  depth  of  yielding  is  the 
limiting  curve  over  most  of  the  contact  loading  space  for  the  given  fi  and  ab.  These  curves  can 
be  generated  quite  easily  from  the  figures  provided  in  the  previous  two  sections  for  many  other 
combinations  of  fi  and  ab)  however  only  two  examples  will  be  displayed  in  this  work.  In  Figure 
5-14  the  optimum  peening  depth,  z/a ,  was  again  found  for  a  contact  load  of  (Po/cry)  =  0.8. 
In  this  example  the  optimum  peening  depth  was  found  to  be;  z/a  ~  0.36.  The  shot  peening 
depth  found  in  this  example  would  be  the  minimum  recommended  because  material  yielding  is 
still  occurring  at  this  depth.  The  material  is  just  at  the  yield  point,  but  it  is  recommended 
that  the  material  be  shot  peened  to  a  slightly  greater  depth  to  ensure  no  material  yielding  and 
hence,  no  mechanical  residual  stress  relieving  will  occur. 
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Figure  5-13:  Optimum  peening  depth  plotted  versus  maximum  contact  pressure  for  fx  —  1.0 
and  at,  =  0.4cry. 


Figure  5-14:  Optimum  peening  depth  plotted  versus  maximum  contact  pressure  for  \x  =  0.8, 
<7b=0.1(7y. 
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Parameter 

Value 

AKo 

3  MPa^n 

R 

0.1 

a 

50  n m 

Y 

1 

Table  5.1:  Parameters  and  respective  values  for  the  threshold  stress  intensity  range  analysis. 

5.3  Stress  Intensity  Range  Threshold  Analysis 

The  maximum  principal  stress  analysis  is  a  conservative  approach.  An  alternative  to  the 
maximum  principal  analysis  is  an  analysis  based  on  the  threshold  stress  intensity  factor  range 
AKq.  At  stress  intensity  factor  ranges  below  the  threshold,  cracks  either  remain  dormant  or 
grow  at  undetectable  rates  [24].  This  could  be  used  in  place  of  the  maximum  principal  stress 
analysis  to  determine  the  optimum  peening  depth.  Using  the  values  listed  in  Table  5.1,  a 
MATLAB  script  file  was  written  to  determine  the  peening  depth  such  that  AKappiieci  <  AKo- 
The  expression  used  for  A Kappiied  is  given  by: 

AKappiied  =  Y(Aa)^na  (5.3) 


where  Y  is  a  geometrical  factor. 

The  analysis  was  performed  using  erj,  =  0.6 ay,  and  p,  =  0.8.  The  material  stresses  in  the 
^-direction  were  used  to  determine  Aa.  Figure  5-15  shows  the  results  of  the  analysis,  where 
the  required  peening  depth  is  plotted  versus  the  maximum  contact  pressure.  In  this  figure  the 
results  of  the  maximum  principal  stress  analysis  for  the  same  loading  conditions  are  also  given 
for  comparison.  As  can  be  seen  in  the  figure,  the  maximum  principal  stress  analysis  is  the 
more  conservative  analysis. 


5.4  Protective  Coatings 

The  protective  coatings  analysis  in  this  work  focuses  on  the  ability  of  the  coating  to  change 
the  coefficient  of  friction  of  the  system.  Again  a  sphere  on  flat  plane  geometry  is  used  in  this 
section.  A  change  in  the  coefficient  of  friction  has  a  direct  impact  on  the  contact  stresses.  The 
contact  stresses  displayed  in  Equations  2.10  and  2.11  are  shown  again  below  in  Equations  5.4 
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Figure  5-15:  Comparison  of  the  required  peening  depths  for  the  maximum  principal  stress 
analysis  and  the  threshold  stress  intensity  range  analysis. 


and  5.5.  The  stresses  at  the  edge  of  contact,  at  the  surface,  are  given  by: 


<rx  =  P0 


3 


.4  +  1/ 
±— ' ’»  • 


x  =  and 


(5.4) 


{2v  —  l)  37 TVfJ, 

3  ~T~ 


X  =  =po. 


(5.5) 


Therefore,  reducing  the  coefficient  of  friction  from  a  value  of  1.2  to  1.0  will  reduce  the  stresses 
from  ax  =  2.160Po  to  ax  =  1.822 P0  and  ay  =  0.291Po  to  ay  =  0.220Po,  assuming  v  —  0.3.  So 
it  can  be  seen  that  changing  the  coefficient  of  friction  has  a  significant  impact  on  the  stresses 
induced  by  contact  and  therefore,  on  the  fretting  fatigue  life  of  the  structure.  Figure  5-5  also 
demonstrates  the  benefit  of  reducing  the  coefficient  of  friction.  It  must  be  remembered  that 
the  equations  shown  above  were  for  global  sliding  conditions.  However,  Chivers  and  Gordelier 
in  [16]  have  shown  that  for  /i  >  0.2,  reducing  fi  will  also  reduce  stress  for  stick-slip  conditions. 

No  experiments  have  been  conducted  as  part  of  this  work  in  which  the  coefficient  of  friction 
was  purposely  altered  to  examine  the  actual  effect  on  fretting  fatigue  life.  However,  Chivers 


and  Gordelier  [41]  performed  a  number  of  experiments  for  various  fretting  fatigue  palliatives 
including  shot  peening  and  protective  coatings  and  their  results  do  support  the  validity  of  the 
discussion  above.  They  examined  two  types  of  non-metallic  coatings.  Both  were  graphite-filled 
resins,  one  of  the  resins  was  a  polyimide  and  the  other  was  an  isocyanate  epoxy.  The  improve¬ 
ments  in  fretting  fatigue  behavior  produced  by  the  coatings  were  quite  significant  and  were 
observed  at  all  stress  levels.  Chi  vers  and  Gordelier  [41], define  a  term  called  the  improvement 
factor,  If,  to  measure  the  success  of  the  various  palliatives  being  investigated.  If  is  given  by: 

If  =  (5.6) 

<J  —  (Tf 

where  crp  is  the  run-out  stress  under  fretting  conditions  with  a  palliative,  cr f  is  the  run-out 
stress  under  fretting  conditions  and  no  palliative  and  cr  is  the  run-out  stress  without  fretting. 
Olivers  and  Gordelier  define  the  run-out  stress  as  the  stress  at  which  failure  does  not  occur 
in  about  108  cycles.  The  improvement  factor  values  were  found  to  be  If  —  0.3  for  graphite 
in  polyimide  resin  and  //  =  0.55  for  graphite  in  epoxy  resin.  For  comparison  purposes  the 
improvement  factors  found  for  shot  peening  were  If  —  0.47  for  low  intensity  peening  and 
If  =  0.43  for  high  intensity  peening.  It  is  interesting  that  the  graphite  in  epoxy  resin  coating 
actually  outperformed  shot  peening.  Olivers  and  Gordelier  also  performed  experiments  using 
a  copper  shim  and  the  combination  of  a  copper  shim  and  graphite  in  epoxy.  The  use  of  the 
copper  shim  was  intended  to  reduce  the  contact  stresses  at  the  surface  of  the  material  due  to 
fretting,  which  from  the  previous  analysis,  again  should  improve  fretting  fatigue  behavior.  The 
improvement  factors  found  were  If  =  0.63  with  copper  shim  and  If  =  0.77  using  both  a  copper 
shim  and  graphite  in  epoxy.  Flat  ended  cylindrical  pads  made  of  ICr-Mo  steel  were  fretted 
against  rectangular  3.5Ni-Cr-Mo-V  specimens  for  the  experiments  conducted  by  Olivers  and 
Gordelier  in  [41].  The  reduction  in  fatigue  life  due  to  fretting  is  directly  related  to  the  stresses 
at  the  contact  surface  and  if  the  stresses  can  be  reduced  the  fatigue  life  will  be  improved. 


Chapter  6 


Conclusions  and  Future  Work 


6,1  Conclusions 

The  key  findings  in  this  work  are  summarized  below. 

•  The  surface  stresses  due  to  contact  can  be  used  to  characterize  the  damaging  effects  of 
fretting  on  the  fatigue  properties  of  a  material.  The  critical  fretting  parameters,  such 
as:  the  normal  contact  load,  the  oscillating  tangential  load  and  the  coefficient  of  friction 
are  all  incorporated  into  the  expressions  for  the  contact  stresses.  It  was  shown  that  for 
a  given  bulk  loading  condition,  the  fretting  fatigue  life  is  reduced  as  the  effective  contact 
stresses  are  increased. 

•  A  fretting  fatigue  knock-down  factor  curve  has  been  produced  for  A1  7075-T6.  It  shows 
the  approximate  relationship  between  the  effective  stress  due  to  contact  and  the  magnitude 
of  the  fretting  fatigue  knock-down  factor,  KF. 

•  It  was  shown  that  the  use  of  a  single  correction  factor  to  modify  the  Basquin  relationship 
for  fretting  fatigue  would  require  extensive  experimental  work.  The  correction  factor  can 
be  expressed  as  a  function  of  the  effective  stress  due  to  contact  for  a  given  bulk  loading 
condition.  However,  this  relationship  changes  as  the  bulk  loading  is  changed.  Therefore, 
families  of  correction  factor  curves  would  need  to  be  developed  for  various  bulk  loadings. 
The  experimental  work  required  may  be  tenable  if  the  correction  factor  was  determined 
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to  be  a  function  of  the  contact  stresses  only  and  not  dependent  on  the  contact  geometry 
for  a  given  material. 

•  A  failure  envelope  curve  has  been  generated  for  A1  7075-T6.  This  curve  shows  the 
approximate  combination  of  bulk  and  contact  loading  that  will  produce  component  failure. 

•  The  damage  tolerant  approach  analysis  has  been  used  to  bound  the  crack  initiation  and 
short  crack  propagation  stages  of  fretting  fatigue. 

•  A  new  method  has  been  proposed  to  determine  the  optimum  peening  depth  for  a  com¬ 
ponent  that  undergoes  fretting  using  the  stress  analysis  provided  by  Hamilton  [14].  The 
method  was  applied  to  a  sphere  on  flat  surface  contact  geometry,  but  is  valid  for  any  con¬ 
tact  geometry  that  has  analytic  expressions  for  the  stress  fields.  The  optimum  peening 
depth  has  been  determined  for  Ti-6A1-4V  for  various  combinations  of  bulk  and  contact 
loadings,  for  a  given  coefficient  of  friction.  This  analysis  was  performed  for  four  different 
values  of  the  coefficient  of  friction,  so  that  the  effect  of  coatings  on  peening  depth  could 
be  analyzed. 

6.2  Future  Work 

More  experimental  work  needs  to  be  performed  to  determine  a  more  accurate  failure  envelope 
curve.  This  could  prove  to  be  a  very  beneficial  design  tool,  especially  if  it  was  shown  to  be 
geometry  independent.  Theoretical  work  is  also  needed  for  this  curve  to  better  understand  why 
the  given  combination  of  contact  and  bulk  stresses  lead  to  component  failure.  This  will  require 
a  thorough  examination  of  the  crack  initiation  stage  of  fretting  fatigue  crack  growth. 

The  shot  peening  analysis  provided  in  this  work  also  needs  experimental  validation  to  de¬ 
termine  the  actual  depth  of  mechanical  stress  relieving  due  to  the  contact  stresses.  If  the 
compressive  residual  stresses  due  to  shot  peening  are  relieved,  the  damage  to  the  surface  by  the 
shot  peening  process  may  cause  a  significant  reduction  in  the  fatigue  strength  of  the  material 
as  was  discussed  earlier. 

More  experimental  work  also  needs  to  be  done  to  determine  the  exact  location  of  crack 
initiation.  The  pictures  displayed  in  the  previous  work  do  not  provide  conclusive  evidence  for 
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the  location  of  crack  initiation.  Interrupted  tests  are  recommended  because  they  could  provide 
not  only  the  location  of  crack  initiation,  but  could  also  provide  information  on  the  number  of 
cycles  required  for  crack  initiation.  More  analysis  is  also  required  for  the  predicted  location  of 
crack  initiation  for  fretted  components.  The  stress  analysis  approach  does  not  always  predict 
the  correct  location  of  crack  initiation. 
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Appendix  A 


Ti-6A1-4V  Yielding  Analysis 
Program 


%Ti-6Al-4V  Yielding  Analysis 
%Input  the  material  properties 
E=115.7e9; 
nu=0.32i; 

R— 12.7e-3; 

Sy=925e6; 

K=3*(l-nu~2)/(2*E); 

Sres=-0.85*Sv 

Saxial=-.l*Sv 

%Where  E: Young’s  Modulus,  nu:Poisson’s  Ratio,  R:Radius  of  fretting  pads, 

%Sy:Yield  Strength,  Sres:Residual  Stress,  SaxiahApplied  bulk  stress 

for  j j = 1 ;  8 

mu=0.5; 

for  k=l:6 

P=0.20; 

Ry=0; 

%  Increase  the  contact  load,  P,  for  each  mu  and  Saxial  until  the  effective 


%  stress  exceeds  Sy. 
while  Ry<l 
a=(K*P*R)~(l/3); 

Po=(3*P)/(2*pi*a~2); 

Bl=Po/a; 

B2=mu*Bl; 

%Input  Hamilton’s  Equations 
y=o; 

z(l)=0.0000001*a; 
x(l)=-1.05*a; 
for  j=l:l 
for  i=l  :1 10 
r=x(i); 

A=(r~2)+(z(j)-2)-(a~2); 

S=((A-2)+(4*(a-2)*(zO)"2)))-(l/2); 

M=((S+A)/2)-(l/2); 

X=((S-A)/2)-(l/2); 

phi=atan(a/M); 

G=(M~2)-(X~2)+(z(j)*M)-(a*N); 

H=(2*M*X)+(a*M)+(z(j)*N); 

%Stresses  due  to  normal  loads 

Sxn=(Po/a)*((l+nu)*z(j)*plu+(l/r''2)*((0-x(i)~2/r''2)*((l-nu)*N*z(j)''2-  . 

((l-2*nu)/3)*(X*S+2*A*N+a"3)-nu*M*z(j)*a)-N*(x(i)"2)-(M*x(i)"2*z(j)*a/S))); 

Syn=(Po/a)*((l+nu)*z(j)*phi-l-(l/r"2)*((x(i)^2/r''2)*((l-nu)*X’*z(j)''2-((l-2*nu)/3) 

*(N*S+2*A*N+a~3)-nu*M*z(j)*a)-N*(2*nu*x(i)''2)-(0/S))); 

Szn=Bl*(-N+(a*z(j)*M)/S); 

TAzxn=Bl*(-z(j)*(x(i)*N/S-(x(i)*z(j)*H/(G~2+H''2)))); 

%Stresses  due  to  tangential  loads 

Sxt=B2*(-x(i)*(nu/4+l)*phi+(a*x(i)*M/r''4)*((3/2-2*x(i)'v2/r''2)*(S*nu-2*A*nu+ 

z0)''2)+(x(i)''2*z(j)''2/S)+(7*nu*r~2/4)-(2*nu*x(i)~2)+r''2)+(x(i)*z(j)*X/r~4)* 


((3/2-2*x(i)~2/r~2)*(-S/6*(l-2*nu)-A/3*(l-2*nu)-.5*(z(j)~2+3*a''2))+ 

(a-2*x(i)-2/S)-(nu*r-2/4)-(7*r-2/4))+(4*a'3*x(i)*z(j)/(3*r^4))* 

(3/2-2*x(i)~2/r~2)*(l-2*nu)); 

Syt=B2*((-3*nu*x(i)*phi/4)+(a*x(i)*M/r''4)*((l/2-2*0/r''2)*(S*nu-2*A*nu+z(j)''2+ 

r~2*nu)+(0*z(j)~2/S)+(3*nu*r~2/4))+(x(i)*z(j)*N/r''4)*((l/2-2*0/r''2)* 

(-S/6*(l-2*nu)-A/3*(l-2*nu)-.5*(z(j)''2+3*a''2))+(a~2*0/S)-(3*nu*r''2/4)-(l*r''2/4))+ 

(4*a''3*x(i)*z(j)/(3*r~4))*(l/2-2*0/r~2)*(l-2*nu)); 

Szt=B2*((z(j)*x(i)*N/(2*r''2)*(l-(r"2+z(j)~2+a''2)/S))); 

TAzxt=B2*((3*z(j)*phi/2)+(a*z(j)*M/r~2)*(l+x(i)~2/r''2-x(i)''2/S)+N/r~2* 

(-.75*(S+2*A)+z(j)/'2-.75*a''2-.25*r''2+z(j)/'2/2*(.5-2*x(i)"2/r''2))); 

%Now  combine  the  stresses 
Sxtot— Sxn+Sxt-{-Sres+Saxial; 

Sytot— Syn+Syt+Sres; 

Sztot=Szn+Szt; 

TAzx  t  ot =TAzxn+TAzxt ; 

%Determine  J2  and  von  Mises  effective  stress 

J2=l  /6*  ( (Sxtot-Svtot) ~  2+(Sytot-Sztot) "  2+  (Sztot-Sxtot) "  2)  +TAzxtot ~  2; 

H2=J2".5; 

Se(i,j)=3-.5*H2; 

x(i+l)=x(i)+.02*a; 

end 

z(j+l)=z(j)+.025*a; 

end 

[Sem,i]=max(Se); 

xl=x/a; 

xl=xr; 

zl=z’; 

zl=zl/a; 

Se=Se; 

%Now  give  x,z  location  of  max  stress 


[Smm ,  zm]  =max  ( Sem) ; 
zmm=z(zm); 
zma— zmm/a; 
xm=i(zm); 
xmm=x(xm); 
xma— xmm/a; 

%Now  find  the  stresses  at  x=a,  and  x=-a  at  the  surface 
a=(K*P*R)~(l/3); 

Po=(3*P)/(2*pi*a^2); 

Bl=Po/a; 

B2=mu*Bl; 

x2=-a; 

r=x2; 

Sxtna=B2*(-pi*x2/2*((nu/4)+l)); 

Sytna==B2*(-3*pi*nu*x2/8); 

Sxnna==Bl*((l-2*nu)*a/3); 

Synna=B  1  *  ( ( l-2*nu)  *  (-a/3) ) ; 

Sxtona=Sxtna+Sxnna+Sres+Saxial; 

Sy  tona= Sy  tna+  Sy nna+  Sres ; 

J  2na=l  /6*  ( (Sxtona-Sy  tona)  ~  2+(Sytona-0)  ~  2+  (O-Sxtona) "  2) +0  ~  2 ; 
H2na=J2na'\5; 

Sena=3~.5*H2na; 

%Find  stresses  at  x=a 

x2=a; 

r=x2; 

Sxtpa=B2*(-pi*x2/2*((nu/4)+l)); 

Sytpa=B2*  (-3*pi*nu*x2/8) ; 

Sxnpa=Bl*((l-2*nu)*a/3); 

Synpa=Bl*((l-2*nu)*(-a/3)); 

SxtoPa=Sxtpa+Sxnpa+Sres+Saxial; 


SytoPa=Sytpa+Synpa+Sres; 

J2pa=l/6*((SxtoPa-SytoPa)~2+(SytoPa-0)~2+(0-SxtoPa)'‘'2)+0~2 

H2pa=J2pa"\5; 

Sepa=3  ~  .5*H2pa; 

%No\v  determine  the  maximum  effective  stress 
Smmm=  [Smm,  Sen  a,  Sep  a] 

[Smmax ,  xmmax]  =max  ( Smmm) ; 

Semax=Smmax; 

Ry=Semax/Sy; 

P=P+.2; 

end 

Pl=P-.2; 

P2(kjj)=Pl; 

mu=mu-hl; 

end 

Saxial= Saxial + ( .  1  *  Sy ) 
end 

mu-mu; 

P2=P2 

%Now  put  P  terms  over  yield  strength, 
for  n=l:8 
for  i=l:6 

a=(K*P2(i,n)*R)~(l/3); 

Pol  (i?n)==(3*P2(i,n))/(2*pi*a~2); 

Ryy(i,n)=Pol(i,n)/Sy; 

end 

end 

Ryy=Ryy 


Appendix  B 


Depth  of  Yielding  Analysis 

%Ti-6Al-4V  Depth  of  Yielding  Analysis 
%Input  the  material  properties 
E=115.7e9; 
nu=0.321; 

R=12.7e-3; 

Sy=925e6; 

K=3*(l-nu"'2)/(2*E); 

Sres=-.85*Sy 

mu=0.3; 

Saxial=-0.1*Sy 
for  jj=l:4 
P=10; 
for  k=l:2 

a=(K*P*R)''(l/3); 

Po=(3*P)/ (2*pi*a~2); 

Bl=Po/a; 

B2=mu*Bl; 

%Input  Hamilton’s  Equations 

y=0; 

z(l)=0.0000001*a; 


x(l)=-1.05*a; 
for  j— 1:27 
for  i— 1:110 
r=x(i); 

A=(i-'2)+(z(j)''2)-(a''2); 

S=((A-2)+(4*(a-2)*(z(j)-2)))-(l/2); 

M=((S+A)/2)-(l/2); 

X=((S-A)/2)-(l/2); 

phi=atan(a/M); 

G=(M-N2)-(X''2)+(z(j)*M)-(a*N); 

H=(2*M*N)+(a*M)+(z(j)*N); 

%STRESSES  DUE  TO  NORMAL  LOADS 

Sxn=(Po/a)*((l+nu)*z(j)*phi+(l/r~2)*((0-x(i)~2/r'‘2)*((l-nu)*N*z(j)''2- 

((l-2*nu)/3)*(N*S+2*A*N+a'‘3)-nu*M*z(j)*a)-N*(x(i)~2)-(M*x(i)''2*z(j)*a/S))); 

Syn^Po/a)*((l+nu)*z(j)*phi+(l/r''2)*((x(i)''2/r~2)*((l-nu)*N*z(j)"2- 

((l-2*nu)/3)*(X*S+2*A*N+a''3)-nu*M*z(j)*a)-N*(2*nu*x(i)''2)-(0/S))); 

Szn=Bl*(-N+(a*z(j)*M)/S); 

TAzxn=Bl*(-z(j)*(x(i)*N/S-(x(i)*z(j)*H/(G~2+H~2)))); 

%Stresses  due  to  tangential  loads 

Sxt— B2*(-x(i)*(nu/4+l)*phi+(a*x(i)*M/r~4)*((3/2-2*x(i)~2/r''2)* 

(S*nu-2*A*nu+z(j)''2)+(x(i)"2*z(j)''2/S)+(7*nu*r~2/4)-(2*nu*x(i)'''2)+r''2)+ 

(x(i)*z(j)*X/r/'4)*((3/2-2*x(i)''2/r''2)*(-S/6*(l-2*nu)-A/3*(l-2*nu)- 

.5*(z(j)''2+3*a~2))+(a~2*x(i)''2/S)-(nu*r~2/4)-(7*r'‘‘2/4))+ 

(4*a''3*x(i)*z(j)/(3*r~4))*(3/2-2*x(i)~2/r~2)*(l-2*nu)); 

Syt=B2*((-3*nu*x(i)*phi/4)+(a*x(i)*M/r^4)*((l/2-2*0/r''2)* 

(S*nu-2*A*nu+z(j),'2+r~2*nu)+(0*z(j)~2/S)+(3*nu*r''2/4))+(x(i)*z(j)*N/r~4)* 

((l/2-2*0/r~2)*(-S/6*(l-2*nu)-A/3*(l-2*nu)-.5*(z(j)~24-3*a~2))+(a~2*0/S)- 

(3*nu*r''2/4)-(l*r''2/4))+(4*a''3*x(i)*z(j)/(3*r''4))*(l/2-2*0/r~2)*(l-2*nu)); 

Szt=B2*((z(j)*x(i)*N/(2*r-'2)*(l-(r''2-fz(j)''2+a''2)/S))); 

TAzxt=B2*((3*z(j)*phi/2)+(a*z(j)*M/r~2)*(l+x(i)''2/r~2-x(i)~2/S)+N/r''2* 

Sfi 


(-.75*(S+2*A)+z(j)~2-.75*a~2-.25*r~2+z(j)~2/2*(.5-2*x(i)~2/r~2))); 

%Now  combine  the  stresses 

Sxtot=Sxn+Sxt+Sres-{-Saxial; 

Sy  tot =Sy  n  4*Sv  ■ t + S  res; 

Sztot=Szn+Szt; 

TAzx  tot  =TAzxn+TAzxt ; 

%Determine  J2  and  von  Mises  effective  stress 

J2=l  /6*((Sxtot-Sytot)"2+(Sytot-Sztot)  "2+(Sztot-Sxtot)  ~2)+TAzxtot~2 
H2=J2'\5; 

Se(ij)=3".5*H2; 

x(i+l)=x(i)+.02*a; 

end 

z(j+l)=z(j)+.02*a; 

end 

[Sem,i]=max(Se); 

xl=x/a; 

xl=xl’; 

zl2=z’; 

zl3=zl2/a; 

Se=Se; 

%Now  give  XjZ  location  of  max  stress 

[Smm,zm]=max(Sem); 

zmm=z(zm); 

zma=zmm/a; 

xm=i(zm); 

xrnm— x(xm); 

xma=xmm/a; 

%Now  determine  when  von  Mises  effective  stress  equals  yield  stress 

Ry=2; 


zma=zma 


zl=zma*a; 
xl(l)=-1.05*a; 
while  Ry>l 
for  i=l:llO 
a=(K*P*R)~(l/3); 

Po=(3*P)/(2*pi*a'"2); 

Bl=Po/a; 

B2=mu*Bl; 

r=xl(i); 

A=(r~2)+(zl~2)-(a~2); 

S=((A-2)+(4*(a-2)*(zl-2)))-(l/2); 

M=((S+A)/2)"(l/2); 

X=((S-A)/2)-(l/2); 

phi=atan(a/M); 

G=(M"2)-(X"2)+(zl*M)-(a*X); 

H=(2*M*X)+(a*M)+(zl*N); 

%Stresses  due  to  normal  loads 

Sxn=(Po/a)*((14-nu)*zl*phi+(l/r~2)*((0-xl(i)~2/r''2)*((l-nu)*X*zl''2- 

((l-2*nu)/3)*(X*S+2*A*N+a''3)-nu*M*zl*a)-N*xl(i)''2-(M*xl(i)"'2*zl*a/S))); 

Syn=(Po/a)*((l+nu)*zl*phi+(l/r~2)*((xl(i)''2/r''2)*((l-nu)*N*zl''2- 

((l-2*nu)/3)*(X*S+2*A*N+a~3)-nu*M*zl*a)-N*(2*nu*xl(i)~2)-(0/S))); 

Szn=Bl*(-N+(a*zl*M)/S); 

TAzxn=Bl*(-zl*(xl(i)*N/S-(xl(i)*zl*H/(G''2+H-'2)))); 

%Stresses  due  to  tangential  loads 

Sxt=B2*(-xl(i)*(nu/4+l)*phi+(a*xl(i)*M/r''4)*((3/2-2*xl(i)~2/r~2)* 

(S*nu-2*A*nu+zl''2)+(xl(i)~2*zl~2/S)+(7*nu*r~2/4)-(2*nu*xl(i)''2)+r~2)+ 

(xl(i)*zl*N/r~4)*((3/2-2*xl(i)~2/r''2)*(-S/6*(l-2*nu)-A/3*(l-2*nu)-.5* 

(zr'2+3*a''2))+(a-'2*xl(i)-'2/S)-(nu*r''2/4)-(7*r''2/4))+(4*a''3*xl(i)*zl/(3*r''4))* 

(3/2-2*xl(i)~2/r~2)*(l-2*nu)); 

Syt=B2*((-3*nu*xl(i)*phi/4)+(a*xl(i)*M/r~4)*((l/2-2*0/r''2)*(S*nu-2*A*nu+zl~2+ 


r~2*nu)+(0*zl~2/S)+(3*nu*r~2/4))+(xl(i)*zl*X/r~4)*((l/2-2*0/r~2)*(-S/6* 

(l-2*nu)-A/3*(l-2*nu)-.5*(zl,'2+3*a''2))+(a''2*0/S)-(3*nu*r''2/4)-(l*r''2/4))+ 

(4*a~3*xl(i)*zl/(3*r~4))*(l/2-2*0/r‘''2)*(l-2*nu)); 

Szt=B2*((zl*xl(i)*N/(2*r''2)*(l-(r''2+zl''2+a''2)/S))); 

TAzxt=B2*((3M*pM/2)+(a*zl*M/r~2)*(l4*l(i)~2/r~2-xl(i)~2/S)+N/r~2* 

(-J5*(S+2*A)+zl~2-J5*a~2-.25*r~2+zl~2/2*(.5-2*xl(i)~2/r~2))); 

%Now  combine  the  stresses 
Sxtot=Sxn+Sxt+Sres+Saxial; 

Sy  tot =Sv  n 4- Sy  t -h  Sres; 

Sztot=Szn+Szt; 

TAzxtot =TAzxn+TAzxt ; 

%Determine  J2  and  von  Mises  equivalent  stress 

J2=l/6*((Sxtot-Sytot)~2+(Sytot-Sztot)~2-f(Sztot-Sxtot)~2)+TAzxtot~2; 

H2=J2~.5; 

Sel(i)=3~.5*H2; 

xl(i-bl)=xl(i)-h02*a; 

end 

Seml=max(Sel); 

Ry=Seml/Sy 

zl=zl4-.01*a; 

end 

zll(k,jj)=(zl/a)-.01; 

P=P+40; 

end 

Saxial=Saxial+.l*Sy 

end 

zl  1— zll 


KQ 


mu=mu; 
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Maximum  Principal  Stress  Analysis 


%Ti-6Al-4V  Maximum  Principal  Stress  Analysis 
%Input  the  material  properties 
E=115.7e9; 
nu=0.321; 

R=12.7e-3; 

Sy=925e6; 

mu=0.30; 

K=3*  (l-nu~  2)  /  (2*E) ; 

Sres=-.85*Sv; 

Saxial=-.l*Sy; 

for  cc=l:8 

P=10; 

for  k=l:ll 

Slmax=l; 

a=(K*P*R)~(l/3); 

Po=(3*P)/(2*pi*a^2); 

Bl=Po/a; 

B2=mu*Bl; 

%Input  Hamilton’s  Equations 


z=0. 0000001  *a; 
x(l)=-1.05*a; 
while  Slmax>0 

%Slmax  is  the  maximum  principal  stress 

for  i=l:110 

r=x(i); 

A=(r-'2)+(z/'2)-(a'2); 

S=((A-2)+(4*(a-2)*(z-2)))-(l/2); 

M=((S+A)/2)~(l/2); 

X=((S-A)/2)-'(l/2); 

phi=atan(a/M); 

G=(M~2)-(X~2)+(z*M)-(a*N); 

H=(2*M*X)+(a*M)+(z*X); 

%STRESSES  DUE  TO  NORMAL  LOADS 

Sxn=(Po/a)*((l+nu)*z*phi+(l/r''2)*((0-x(i)~2/r''2)*((l-nu)*N*z''2- 
((l-2*nu)/3)*(X*S+2*A*N+a''3)-nu*M*z*a)-X*(x(i)~2)-(M*x(i)''2*z*a/S))); 
Syn=(Po/a)*((14-nu)*z*phi+(l/r'‘2)*((x(i)'‘2/r''2)*((l-nu)*N*z''2- 
((l-2*nu)/3)*(X*S+2*A*N+a"3)-nu*M*z*a)-X*(2*nu*x(i)"2)-(0/S))); 
Szn=Bl*(-N+(a*z*M) /S) ; 

TAzxn=Bl*(-z*(x(i)*X/S-(x(i)*z*H/(G''2+H''2)))); 

%Stresses  due  to  tangential  loads 

Sxt=B2*(-x(i)*(nu/4+l)*phi+(a*x(i)*M/r''4)*((3/2-2*x(i)''2/r''2)* 

(S*nu-2*A*nu+z''2)+(x(i)~2*z''2/S)+(7*nu*r''2/4)-(2*nu*x(i)~2)+r-'2)+ 

(x(i)*z*N/r''4)*((3/2-2*x(i)~2/r~2)*(-S/6*(l-2*nu)-A/3*(l-2*nu)- 

.5*(z~2+3*a~2))+(a~2*x(i)''2/S)-(nu*r~2/4)-(7*r~2/4))+(4*a''3*x(i)* 

z/(3*r~4))*(3/2-2*x(i)~2/r~2)*(l-2*nu)); 

Syt=B2*((-3*nu*x(i)*phi/4)+(a*x(i)*M/r'’4)*((l/2-2*0/r''2)* 

(S*nu-2*A*nu+z~2-l-r''2*nu)+(0*z''2/S)+(3*nu*r''2/4))+(x(i)*z*N/r~4)* 

((l/2-2*0/r''2)*(-S/6*(l-2*nu)-A/3*(l-2*nu)-.5*(z~2+3*a''2))+ 

(a''2*0/S)-(3*nu*r''2/4)-(l*r''2/4))+(4*a''3*x(i)*z/(3*r''4))* 
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(l/2-2*0/r~2)*(l-2*nu)); 

Szt=B2*((z*x(i)*X/(2*r"2)*(l-(r's2+z"2+a''2)/S))); 
TAzxt=B2*((3*z*phi/2)+(a*z*M/r"2)*(l+x(i)^2/r-2-x(i)^2/S)+ 
X/r''2*(-.75*(S+2*A)+z~2-.75*a''2-.25*r~2+z~2/2*(.5-2*x(i)''2/r~2))) 
%Now  combine  the  stresses 
Sx  tot = Sxn + Sx  t + Sres+ Saxial ; 

Sytot=Syn4-Syt+Sres; 

Sztot=Szn+Szt; 

TAzx  t  ot =TAzxn-f  TAzxt ; 

Sxtotv(i)=Sxtot; 

Sztotv(i)=Sztot; 

TAzxtot  v  (i)  =TAzxtot ; 

%Determine  the  max  principal  stress  for  each  point 

cxz=(Sxtot+Sztot)/2; 

rxz = ( (Sxtot-cxz)  ~  2+TAzxtot  ~2)~.5; 

Slxz(i)=cxz-frxz; 

x(i+l)=x(i)+.02*a; 

end 

%Find  the  maximum  principal  stress  of  the  vector  of  principal  stresses 
Slxzm=max(Slxz); 

Sxmax=max(Sxtotv); 

Szmax=max  (Sztot  v) ; 

TAzxmax=max(TAzxtotv) ; 

Slmax=Slxzm 

z=z+.01*a; 

zaa=z/a 

end 

z2(k,cc)=z-.01*a; 

a=a; 

z3(k,cc)=(l/a)*z2(k,cc); 
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P=P+40 

end 

Saxial— Saxial + .  1  *  Sy 

end 

z2=z2 

z3=z3 

mu=mu; 


os 


Appendix  D 


Hamilton’s  Expressions 


D.l  Normal  Load  Stresses 


O'  x  — 


3  P 
2tx  a,3 


(1  +v)z<!>  + 


1  f  [(1  _  „)  Nz2  -  (NS  +  2AN  +  a3)  -  uMza] 

■A  J - / 


Mx^za 


-N  (x2  +  2 uy2)  -  Mg 


(Tv  = 


3  P 
2na3 


1  i  [(i  -p)Nz*-l=giL  [NS  +  2AN  +  a3)  -vMza) 

(1  +  v)zd>+-pi\  -N(y2  +  2vx2  )-Mgsz 


3  P 


Oy  ~ 


azM 


z  27t  a3 

At  y  =  0,  the  only  non- zero  shear  stress  is  rzx . 

3P 


T rx  — 


27ra3 


LJ^ 

v 

1  ^ 

G2  +  H 2  J 

where 


A  =  r2  4-  z2  -  a2;  5  =  (A2  +  4a2£2)1/2;  r2  —  x2  +  y 2  and 
m=(S±4)1/2;  7v=(^)1/2;  0  =  tan-1(^)  and 


G  =  M2  -  N2  +  zM  -  aN ;  H  =  2A/JV  +  aM  +  ziV 

On  the  axis  where  r  =  0,  the  formulae  above  are  indeterminate  and  must  be  re-worked. 
Hamilton  performed  this  and  the  expressions  are  now  given  by: 

a3 

(1  +  u)  (z  tan  1  (a/z)  -  a)  +  2^q2  +  ^ 


az  = 


3  P 
27r  a3 


a 1  +  zl 


and  t xy  —  ^~yz  —  t zx  —  0* 


D.2  Tangential  Load  Stresses 


(J  X  — 


3  Q 
27ra3 


-*  (f  +  1)  0  +  {  (|  -  *£)  (Sv  -  2Av  +  z2)  +  ^  -  2ux2  +  r2  } 

+^{(§-^)  [-|(l-2^)-#(l-2z,)-i(z2  +  3a2)]+^-i2f±2d} 

(§-*£)  (1-2./) 

_3^  +  {  (i  -  [(5!/  -  2Av  +  r2v)  +  z2)  +  } 

+igfi  (i  -  ¥)  a " 


Again  the  only  non-zero  shear  stress  at  y  =  0  is  rzx. 

*  Note  the  equation  for  ax  in  the  original  paper  by  Hamilton  contained  an  error.  The  paper 
contained  the  number  3  in  the  expression  marked  above  where  it  should  have  contained  an  S. 
The  units  do  not  work  if  this  change  is  not  made. 


3  Q 
27ra3 


3 zcj)  azM  f  x2  x2>\  N 

^  +  —  {'  +  ?-s}  +  3 


-f  {S  +  2A)  +  z2- 


3  a2+r 
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